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PREFACE 


The  purpose  of  this  study  was  to  examine  the  feasi¬ 
bility  of  using  a  laser  radar  or  lidar  to  track  aircraft 
from  space.  As  a  first  step  to  a  complete  analysis,  this 
thesis  investigated  the  laser  energy  propagation  as  defined 
by  the  lidar  range  equation.  The  results,  generated  from  a 
computer  model,  indicate  that  aircraft  could  potentially  be 
tracked  from  a  satellite  in  low-earth  orbit.  Further  study 
is  necessary  to  analyze  more  atmospheric  conditions  and 
specific  hardware  configurations. 

I  had  a  great  deal  of  help  from  a  number  of  people.  I 
wish  to  thank  Mr.  Frank  Jenks,  Mr.  Ron  Rodney,  Major  Bob 
Hughes,  and  Mrs.  Kris  Larsen  for  their  time,  insight  and 
patience  in  getting  FASCOOE  to  operate.  Also,  I  wish  to 
thank  Captain  Steve  Spence  for  his  superb  insight  ana  help 
on  making  the  computer  model  "user-friendly”.  Also,  many 
thanks  go  out  to  Captain  Don  Holland  and  Dr.  Paul  McManamon 
for  their  assistance  and  insights  which  added  immeasurably 
to  this  thesis.  I  am  indebted  to  and  wish  to  thank  Major 
T.S.  Kelso  (my  reader)  and  Major  David  Stone  (my  advisor) 
for  their  assistance,  corrections,  suggestions  and  patience. 
Above  all,  I  wish  to  thank  my  wife  Pattie  and  my  children, 
John,  Becky,  and  Amanda  for  their  understanding  and 
sacrifice  during  the  many  months  of  work  on  this  thesis. 
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ABSTRACT 


/ 

y 

This  study  is  a  general  analysis  of  the  feasibility  of 
tracking  aircraft  from  space  using  a  laser  radar  or  lidar. 
The  basis  for  the  analysis  is  the  lidar  range  equation. 
Specific  hardware  is  not  discussed. 

A  computer  program  was  written  to  model  the  different 
parameters  involved  with  the  analysis.  Atmospheric  atten¬ 
uation  equations,  which  were  developed  from  FASCODE  data, 
define  the  atmospheric  losses  for  two  type  of  atmospheres 
(clear  and  haze).  Other  lidar  range  equation  variables  are 
either  computed  from  input  data  or  are  directly  inserted 
into  the  program  for  analysis.  Four  different  types  of 
lasers  were  investigated.  The  output  data  from  the  model 
was  plotted  to  show  the  relationship  of  various  input 
parameters  with  the  power  available  at  the  receiver.  The 
results  indicate  that  lidar  systems  carried  on  low-earth 
orbiting  satellites  may  have  the  potential  to  track  air¬ 
craft.  Further  analysis  is  required  to  develop  complete 
feasibility  windows  over  a  wider  range  of  conditions. 
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ANALYSIS  OF  SPACE-BASED 
LIDAR  FOR  AIRCRAFT  TRACKING 


I  -  Introduction 


Tgpic 

This  thesis  is  a  study  of  the  feasibility  of  using  laser 
radar  or  lidar  to  track  aircraft  from  space.  For  this 
study,  the  term  lidar  will  refer  to  an  active  sensor  system 
in  which  a  laser  is  used  to  illuminate  a  target  and  the 
reflection  of  the  laser  energy  from  that  target  is  collected 
and  analyzed.  The  analytical  portion  of  this  study  is 
centered  around  computer  modeling  of  several  key  design 
parameters  of  a  lidar  system. 

Background  &  Justification 

The  news  in  the  United  States  during  the  1980’s  has  been 
replete  with  stories  about  aircraft  accidents.  Many  of 
these  accidents  occurred  due  to  the  lack  of  proper  or 
sufficient  tracking  and  monitoring  of  the  aircraft.  As  the 
number  of  aircraft  flights  around  the  globe  continues  to 
increase,  so  will  the  need  for  tracking. 

Agencies  interested  in  tracking  aircraft  in  the  United 
States  include  the  Department  of  Defense  and  the  Federal 
Aviation  Administration.  Whether  in  the  Persian  Gulf  or  on 
the  Atlantic  coast,  the  Department  of  Defense  is  concerned 
with  maintaining  a  watchful  eye  on  aircraft  that  may 


threaten  the  United  States  or  U.S.  interests  and  forces 
abroad.  The  Federal  Aviation  Administration  meanwhile, 
monitors  the  skies  of  the  United  States  to  provide  for  the 
safety  of  air  travelers  and  the  general  populous  on  the 
ground.  This  monitoring  function  is  heavily  dependent  on 
precise  knowledge  of  aircraft  positions. 

Internationally,  airline  travel  over  the  broad  expanses 
of  ocean  poses  challenges  to  aircraft  traffic  control.  The 
lack  of  precise  position  data  by  an  international  airliner 
could  result  in  a  catastrophe  such  as  an  errant  aircraft 
wandering  into  hostile  airspace. 

Radar  has  been  the  classical  answer  to  the  tracking 
problem  in  the  past.  However,  radar  has  several  limi¬ 
tations.  Depending  on  height  of  the  antennas  and  the  local 
terrain,  ground-based  radar  can  only  detect  and  track 
aircraft  in  the  region  of  the  radar  site.  Low  flying 
aircraft  are  more  difficult  to  detect  due  to  terrain  and 
1 ine-of-sight  limitations.  Ground-based  radar  systems  are 
therefore  very  dependent  on  geography,  both  physical  and 
pol i t i cal . 

As  man’s  quest  for  space  and  his  use  of  satellites 
grows,  new  ideas  for  space-based  sensors  have  arisen.  One 
such  idea  involves  the  use  of  passive  infrared  detection 
technology  (15:1).  Another  space-borne  sensor  concept  also 
found  in  the  literature  is  a  space-based  radar  (5:47). 
Space-based  sensors  have  the  significant  advantage  of  un- 
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precedented  geographical  access,  with  some  orbits  providing 
worldwide  coverage. 

The  function  of  tracking  aircraft  from  space  poses  a 
multitude  of  problems.  The  tracking  system  must  be  able  to 
first  detect  the  aircraft  and  then  to  acquire  data  on  that 
aircraft  over  a  period  of  time  so  as  to  establish  its 
position  and  heading. 

As  one  analysis  showed,  passive  IR  detection  is  limited 
in  position  accuracy  on  the  order  of  1.67  km  of  resolution, 
and  therefore  subject  to  spatial  ambiguities  (15:50). 
Conventional  radar,  operating  in  the  microwave  frequency 
range,  has  a  broad  beamwidth  (dependent  on  antenna  size)  and 
is  therefore  also  susceptible  to  spatial  ambiguities  at  long 
ranges  (8).  To  provide  accurate  position  data  at  long 
ranges,  such  as  from  space,  one  would  like  to  have  a  narrow 
beam  system,  such  as  lidar. 

Scope 

This  thesis  is  limited  to  the  study  of  lidar  to  track 
aircraft  from  space.  Although  other  tracking  systems  and 
techniques  may  be  feasible,  this  analysis  investigates  only 
the  1 i dar  option . 

Further,  this  analysis  is  general  in  nature.  Specific 
details  on  laser  physics,  signal  processing  and  atmospheric 
attenuation  are  omitted.  Detail  in  these  areas  is  the 
subject  of  books  and  would  drive  this  analysis  into 
unreasonable  length  and  complexity.  Therefore,  broad  but 
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realistic  categorizations  of  the  laser,  the  processing  and 
the  atmosphere  are  used. 

Finally,  this  thesis  is  not  an  attempt  to  develop 
hardware  nor  new  technologies.  Technological  capabilities 
in  hardware  are  assumed  and/or  speculated.  The  model  is 
founded  on  general  orbital,  optical  and  lidar  principles  and 
should  remain  useful  even  as  the  technology  advances. 


II .  Literature  Review 


THE  LIDAR  SYSTEM 

All  lidar  systems  can  broken  down  into  three  major 
components:  the  laser,  the  detector  and  the  processor.  When 
designing  a  lidar  system,  the  choice  of  specifics  for  these 
three  major  subcomponents  is  a  function  of  the  mission  and 
tasking  of  the  lidar  system  (i.e.,  what  it  is  supposed  to 
detect).  To  track  aircraft,  the  laser  must  operate  at  a 
wavelength  that  can  be  propagated  through  the  atmosohere. 
Once  a  suitable  wavelength  is  chosen,  a  detector  and  aopro- 
priate  processing  technique  must  be  selected  to  obtain  the 
required  information  from  the  signal  return. 

When  choosing  the  laser  for  a  lidar  system,  two 
important  characteristics  must  be  considered:  1  )  the 
atmospheric  attenuation  at  the  operating  wavelength  and  2; 
the  potential  eye  hazard  associated  with  lasers  operating  at 
that  wavelength. 

Atmospheric  attenuation  is  a  function  of  the  aerosols 
in  the  optical  path  as  well  as  the  type  and  density  of  the 
gases  encountered  along  the  path.  Certain  gases  absorb  at 
certain  wavelengths  while  they  are  relatively  transparent  at 
other  wavelengths.  Likewise,  certain  aerosols  absorb  or 
scatter  stronger  at  some  -wavelengths  than  others.  In  combi¬ 
nation,  the  gases  and  aerosols  in  the  atmosphere  create 
spectral  pockets  of  high  attenuation  and  spectral  windows  of 


high  transmittance.  To  maximize  the  range  of  a  lidar 
system,  a  designer  would  chose  a  laser  that  operates  at  a 
wavelength  that  corresponds  to  one  of  the  atmospheric 
windows.  This  study  looks  at  four  such  laser  wavelengths, 
all  of  which  happen  to  operate  in  an  atmospheric  window  in 
the  infrared  portion  of  the  spectrum. 

The  eye-safe  issue  for  laser  operation  is  again  depend¬ 
ent  on  the  wavelength  of  operation.  Sliney  and  Wolbarsht, 
adapting  from  the  American  Conference  of  Governmental 
Industrial  Hygienists,  break  the  IR  laser  wavelengths  into 
three  categories  (25:261,262,267).  The  first  category,  700 
nm  to  1049  nm,  has  the  most  restrictive  exposure  limits 
while  the  1.4  Mm  to  10^  Mm  band  has  the  least  restrictive 
limits  (25:262).  Essentially,  the  closer  the  wavelength  is 
to  the  visible  wavelengths,  the  more  restrictive  the 
exposure  limits  become. 

Dr.  Paul  McManamon  from  the  Air  Force  Wright  Aero¬ 
nautical  Laboratories,  has  written  on  two  type  of  lasers  for 
potential  use  in  space,  CO^  and  Neodymium  (17:3).  Both 
operate  in  an  atmospheric  window,  however  Neodymium  could  be 
an  eye  hazard  with  an  operating  wavelength  around  1.06  Mm. 
CO^,  operating  at  10.6  Mm,  has  an  exposure  limit  rating  more 
than  1000  times  higher  than  neodymium  (25:262), 

In  light  of  the  eye  hazards  of  neodymium.  Dr.  McManamon 

suggested  two  other  lasers  with  longer  wavelengths,  for 

1  ^  16 

investigation:  Holmium,  operating  at  2.1  Mm,  and  C  '0,  , 
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13  16 

operating  at  11.1  (18).  The  designation  C  0^  refers 

to  carbon  dioxide  (CO^)  distinguished  by  the  use  of 
carbon-13  instead  of  carbon-12  in  the  molecule.  Since  the 
attenuation  from  atmospheric  gases  is  a  function  of  the 
density  of  the  absorbing  gas,  one  wishes  to  avoid  operating 
the  laser  at  a  wavelength  that  corresponds  to-  the  absorption 
bands  of  any  of  the  prevalent  atmospheric  gases.  The 
atmosphere  contains  very  little  carbon-13  based  carbon 
dioxide.  Hence,  atmospheric  attenuation  is  lower  when  using 
a  carbon- 13  based  CO^  laser  than  a  carbon- 12  based  CO^ 
laser. 

The  other  major  portions  of  a  lidar  system  are  the 
detector  and  the  associated  processing  subsection.  The 
actual  choice  of  detector  material  is  wavelength  dependent 
and  often  effective  over  only  a  small  bandwidth.  For 
example,  Kane,  Zhou  and  Byer  recommend  a  silicon  detector 
for  their  Nd:YAG  lidar  system  (13:2481).  Silicon  has  a 
cutoff  wavelength  at  1.1  pm  (4:191).  Therefore,  a  silicon 
detector  would  be  sufficient  for  Nd:YAG  but  not  for  CO,. 

The  processing  involved  With  lidar  systems  is  another 
factor  that  would  be  an  integral  design  consideration  and  a 
function  of  the  mission  of  the  lidar  system.  Since  the 
system  under  study  here  would  be  required  to  track  low 
altitude  aircraft,  ground  returns  must  be  differentiated 
from  targets.  The  targets  could  be  differentiated  from  the 
ground  by  their  doppler  signature.  An  aircraft,  unless 
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traveling  with  no  relative  motion  with  respect  to  the 
ground,  would  impart  a  doppler  shift  on  the  laser  energy 
which  is  different  than  the  ground’s  doppler  shift.  This 
difference,  if  large  enough,  could  then  be  detected. 

One  way  to  accomplish  the  doppler  detection  is  by 
heterodyning.  Heterodyning  is  a  procedure  in  which  a  signal 
is  mixed  with  a  stable  frequency  source  and  the  difference 
between  the  two  frequency  sources  is  extracted.  Fluckiger, 
Keyes  and  Shapiro  describe  this  process  as  occurring  on  the 
photodetector  itself  (10:318). 

Kane,  Zhou  and  Byer  propose  heterodyne  detection  in 
their  design  of  a  doppler  wind  velocity  lidar.  (1.1:2481) 
Their  implementation  of  heterodyning  references  the 
frequency  shifting  of  a  Nd:YAG  signal  by  "acoustoopt i cal  or 
electrooptic  modulation."  03:2480)  They  also  state: 

Pulsed  reference  beam  coherent  lidar  systems  put  two 
principal  requirements  on  the  laser  technology:  first, 
the  'ow-power  oscillator  must  be  stabilized  so  that  the 
frequency  drift  during  the  pulse  round  trip  time 
contributes  only  a  small  fraction  of  the  allowable 
frequency  error;  second,  a  high-energy  pulse  must  be 
transmitted  at  a  known  offset  frequency  from  the 
oscillator  to  well  within  the  allowable  frequency 
error .  (13: 2480 ) 

They  further  state,  with  two  references  listed,  that  such 
frequency  stability  is  achievable  in  Nd:YAG  lasers  operating 
at  low  powers  (13:2480). 

Another  method  of  extracting  doppler  information  from  a 
returned  signal  is  with  autodyne  detection,  as  proposed  by 
Fluckiger,  Keyes  and  Shapiro  (10:318). 
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Autodyne  detection  employs  no  local  oscillator  beam. 

It  is  a  direct  detection  scheme  involving  self-beating 
between  the  various  frequency  components  of  the 
received  signal  beam.  (10:318) 

Such  a  technique  obviates  the  necessity  for  the  extreme 
pulse-to-pulse  frequency  stability  as  discussed  by  Kane, 

Zhou  and  Byer  above  (10:318).  This  specific  point  is  very 
important  with  space-based  lidar  systems  since  the  range  to 
target  is  long,  resulting  in  a  longer  time  base  in  which 
pulse-to-pulse  frequency  stability  is  required.  FlucKiger, 
Keyes  and  Shapiro  further  explain  that  autodyne  detection 
can  be  used  for  doppler  processing,  can  be  as  sensitive  as 
heterodyne  detection"  (10:325)  and  "is  far  less  sensitive  to 
receiver  optics  quality"  (10:325)  than  heterodyne  detection 
( 10:325 ) . 

THE  ENVIRONMENT 

In  addition  to  the  components  of  the  lidar  system,  a 
designer  would  need  to  consider  the  sources  of  attenuation 
that  the  system  would  be  likely  to  encounter.  For  this 
thesis,  the  atmosphere  poses  the  largest  single  source  of 
attenuation.  The  effects  of  attenuation  can  be  mathemat¬ 
ically  modeled  with  the  lidar  range  equation.  The  partic¬ 
ular  version  of  the  lidar  range  equation  used  for  this  study 
is  taken  from  Dr.  McManamon  (17:1): 


where 
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P^  =  power  received  on  the  detector 

P^  =  power  transmitted 

=  transmitter  gain 

=  target  cross-section 

R  =  range  from  the  lidar  to  the  target 

=  area  of  the  receiving  optics  and 

r  =  total  (one-way)  attenuation  from  all  sources 

(17:1) 

Several  atmospheric  attenuation  models  are  available 
for  calculating  the  atmospheric  attenuation.  Two  atmos¬ 
pheric  attenuation  models  were  investigated  for  this  tnesis. 
LOWTRAN  and  FASCOOE.  Both  software  packages  are  products  of 
the  Air  Force  Geophysics  Laboratories.  FASC0D2  is  the 
particular  version  of  FASCOOE  that  was  available  at  the  time 
of  this  thesis  (20)  and  the  terms  are  used  interchangeably 
throughout  this  document. 


The  LOWTRAN  code  calculates  atmospheric  |:ransmit- 
tance  and  radiancy,  averaged  over  20-cm'  intervals 
in  steps  of  5  cm'  from  350  to  40,000  cm  tO.25  to 
,  28.5  fJm).  The  code  uses  a  single-parameter  band 

model  for  molecular  absorption,  and  includes  the 
effects  of  continuum  absorption,  molecular  scatter¬ 
ing,  and  aerosol  extinction.  (14:10) 


FASC0D2  on  the  contrary,  "is  a  model  and  computer  coae  for 
the  accelerated  line  by  line  calculation  of  spectral 
transmittance  and  radiance  for  atmospheric  problems... ' . 
(6:1)  In  other  words,  FASCOOE  is  narrowband.  Since  a  laser 
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operates  at  such  a  narrow  linewidth  or  narrowband.  FASCODE 
appears  to  be  more  appropriate  than  LOWTRAN  for  laser 
attenuation  modeling  (12; 16). 

To  illustrate  this  point,  consider  the  atmospheric 
attenuation  for  a  CO^  laser.  The  laser  energy  would  be 
absorbed  by  the  atmospheric  CO,  at  a  very  specific  wave¬ 
length.  This  specific  absorption  is  modeled  in  FASCODE  in 
its  line-by-line  calculations.  In  LOWTRAN  however,  the 
narrow  band  laser  energy  attenuation  would  be  approximated 
by  the  broadband  average  over  a  wide  bandwidth.  Because  of 
LOWTRAN’s  bandwidth  mismatching  and  the  corresponding 
errors,  FASCODE  was  deemed  more  appropriate  for  laser 
attenuation  modeling. 


Ill . 


The  Model 


The  main  objective  of  this  study  was  to  determine  the 
general  feasibility  of  using  a  lidar  system  on  a  satellite 
to  track  aircraft.  The  prime  tool  in  accomplishing  this 

p 

task  was  a  computer  model.  This  chapter  explains  the 
computer  model  and  many  assumptions  behind  the  model.  The 
code  was  written  in  Microsoft^  QuickBASIC  on  an  IBM.  PC/XT 
compatible  machine.  The  name  given  to  the  model  is  the 
Optical  Aircraft  Tracking  System  (OATS)  model.  The  actual 
computer  code  for  the  model  is  listed  in  Appendix  A. 

The  OATS  model  was  written  to  be  a  flexible  engineering 
and  research  tool  to  investigate  various  design  alterna¬ 
tives,  The  model  does  not  include  all  alternatives  avail¬ 
able  to  a  designer  nor  does  it  delve  into  the  specifics  of 
individual  subsystem  operations.  For  example,  the  laser 
subsystem  is  modeled  with  four  inputs;  1)  the  laser  pulse 
power,  2)  the  wavelength  of  operation  3)  the  laser  pulse 
width  and  4)  the  beam  divergence  angle  for  the  post-optics 
beam.  Specific  laser  mechanics  such  as  the  lasing  medium 
and  the  pulsing  mechanism  are  not  discussed.  A  similar 
level  of  generalization  was  used  for  the  other  subsystem 
models  as  well. 

Equation  1,  the  lidar  range  equation,  is  the  core  of 
the  model.  The  variables  for  the  equation  are  either  given 
as  inputs  or  calculated  from  input  data.  Subprograms  or 
functions  are  used  for  much  of  the  variable  manipulation  and 


generation.  Output  data  is  selectable  from  a  variety  of 
options  available  to  the  modeler.  For  example,  the  modeler 
may  need  the  available  target  track  time  per  satellite 
orbit.  This  is  based  on  the  zenith  angle  at  which  a  target 
is  first  detected.  To  calculate  an  initial  detection  zenith 
angle,  the  model  uses  multiple  iterations  of  the  lidar  range 
equation  to  compute  certain  parameters  and  match  them  witn 
the  initial  detection  criterion.  Once  this  angle  is  found, 
a  separate  function  computes  the  track  time  available  during 
a  single  pass  of  the  satellite. 

The  OATS  model  can  be  broken  into  seven  major  subsets; 

1)  Orbital  Mechanics 

2)  The  Ranging  Algorithm 

3)  The  Transmittance  Algorithm 

4)  The  Target  Lidar  Cross  Section  Algorithm 

5)  The  Lidar  Range  Equation 

6)  The  Inputs 

7)  The  Output  Data  and  Calculations 

Each  of  the  seven  blocks  involve  several  major  assumptions 
and  core  concepts.  These  assumptions  and  concepts  will  be 
discussed  in  the  following  sub-chapters.  Appendix  B 
contains  the  master  list  of  all  variables,  an  explanation  of 
each  variable  and  the  name  given  to  that  variable  in  the 


code . 


Because  the  lidar  system  under  discussion  will  be 
carried  on  a  satellite,  orbital  mechanics  plays  an  important 
role  in  the  OATS  model.  The  reference  for  the  analysis  of 
the  geometry  will  be  the  target.  The  zenith  angle  (from 
the  target  to  the  satellite)  will  be  used  to  vary  the 
location  of  the  satellite.  To  simplify  the  problem  ,  only 
circular  orbits  are  discussed  and  explored.  Elliptical, 
parabolic,  hyperbolic,  or  ballistic  trajectories  are  also 
possible  but  will  not  be  addressed  in  this  analysis. 

Another  simplification  in  the  orbital  mechanics  of  this 
modeling  effort  is  that  the  earth  is  a  perfect  sphere. 
Actually,  the  earth’s  equatorial  radius  is  6378.145  km  and 
the  polar  radius  is  6356.785  km  (3:94).  However,  in  lieu  of 
a  more  complex  model  of  the  earth,  the  OATS  model  uses  a 
spherical  approximation.  The  OATS  model  concentrates  on 
midlatitude  targets  only  and  therefore  uses  the  Midlatitude 
Winter  default  value  found  in  FASCODE  of  6371.230  km  (13:5) 
for  the  earth’s  radius.  This  value  is  assumed  to  be  the 
constant  radius  of  the  earth  for  all  computations.  The 
errors  induced  by  this  approximation  are  very  small,  causing 
a  variation  in  the  signal  power  received  at  the  satellite  of 
less  than  one  percent. 

To  define  the  length  of  time  that  the  satellite  taKes 
to  complete  one  orbit  of  the  earth  (the  orbital  period),  the 
speed  of  the  satellite  must  first  be  known.  The  speed  or 


velocity  of  a  satellite  is  a  function  of  its  altitude.  The 
scalar  representation  of  this  relationship  (taken  from  Bate, 
Mueller  and  White)  is  given  in  Equation  2  as  follows: 


where  n  -  3.986012  x  10® 

sec 

r  =  radius  of  the  orbit  and 
V  •  velocity  of  the  satellite  in  km/sec  (3:34,429). 

S 

The  radius  of  orbit  is  the  distance  from  the  center  of  mass 
of  the  earth  to  the  center  of  mass  of  the  satellite  ana 
approximated  here  to  be  the  sum  of  the  radius  of  the  earth 
and  the  satellite  altitude.  Therefore, 


r  =  R  +  h  i  j ) 

where  R  is  the  radius  of  the  earth  (6371.23  km)  and  h  is 

e 

the  altitude  of  the  satellite  from  the  earth's  surface. 

Using  Equation  2,  one  can  readily  calculate  the  orbital 
period  of  a  satellite.  The  circular  distance  in  one  orbit 
is  simply  the  circumference  of  the  orbit  or  2Tr(  r  ).  The 
orbital  period  therefore,  is  the  circumference  of  the  orbit 


2nr 

^  V 

8 


2nr 


(4) 


/  ' 

2nr  /  r  2^ 

P  =  -  where  -  =  K1  ( 5 ) 


such  that  K\  -  0.009952  when  r  is  in  km  and  P  is  in 
seconds.  Therefore,  the  orbital  period  in  seconds  is 
expressed  as 


P  =  K1 


!  6  I 


From  this  relation,  one  can  easily  calculate  the  angular 
revolution  rate  of  the  satellite.  The  rate.  in  radians 
per  second  is  simply  expressed  as: 


w 


2TT 

P 


or  y  = 


631 .349 


.  6 


(  7  ) 


Given  a  zenith  angle  6  and  the  corresponding  geocentric 

angle  defined  in  Figure  1,  the  geocentric  angular  rate 

d<l>  , 

or  - ■  is  equivalent  to  y. 
dt 

The  geometry  in  Figure  1  yields  a  triangle  for  whicn 
Equation  8  applies  (9:8). 
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FIGURE  I  ORBITAL  GEOMETRY 


R  d  sin(C  ) 


(8) 


2 


h 


r 


From  Figure  1  and  Equation  8,  one  can  calculate  the  transit 
time  from  a  given  zenith  angle  inbound  to  the  same  zenith 
angle  outbound  (i.e.  symmetrical  about  the  normal  to  the 
earth’s  surface)  for  a  circular  orbiting  satellite  at  any 
given  altitude  (assuming  that  the  satellite  passes  through 
the  normal  vector).  In  Figure  1,  d  is  the  distance  from  the 
satellite  to  a  point  or  target  on  earth  while  R  is  the 
radius  of  the  earth  (assumes  that  the  target  is  on  the 
earth’s  surface).  Also,  the  angle  8  defines  the  angle  C 
since  (C  +  8)  must  equal  n.  Therefore,  Equation  8  can  be 
simplified  as ; 


I  8  ) 


(10) 


However,  from  the  given  geometry,  h  is  also  defined  as 

C 


hr=R  dsin(C) 

C  « 

or 

R  d 

h  =  — ^ —  sin  (C) 
c  r 


h  =  R  sin  (^) 

c  e 


(11; 


which  yields: 


R  s  i  n  ( t )  =  — - —  sin  (  C  ) 
e  r 


(12; 
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l  1  4  ) 

Given  the  geocentric  angle  1>  which  is  defined  by  d,  r 
and  0  (where  0  =  n  -  c,  and  d,  r,  and  0  are  entering  argu¬ 
ments),  the  rate  of  change  of  <t>  is  simply  the  satellite’s 
revolution  rate  o.  Given  the  two  dimensional  case  (when  the 
satellite  passes  directly  overhead  of  the  target)  and  given 
that  0^  is  the  zenith  angle  where  tracking  begins,  the  time 
from  0^  to  when  the  satellite  is  directly  overhead  (0=0) 
is  simply  Therefore,  if  tracking  is  symmetric  aoouc 

the  0  =  0  point,  then  the  time  of  potential  track  >p  is 
Note  that  this  assumes  that  the  track  will  not  be  lost 
through  the  doppler  null  (directly  overhead)  or  for  any 
other  reason  such  as  fluctuating  Lidar  Cross  Section  ( lC3  )  . 
which  is  a  function  of  aspect  angle. 

The  Ranging  Algorithm 

To  calculate  range  one  must  first  establish  a  refer¬ 
ence.  In  Figure  2.  the  origin  of  the  three  axes  is  the 
location  of  the  target  and  therefore  the  central  reference 
point.  The  target  here  is  assumed  to  be  at  low  level  i  10 
meters  from  the  surface  for  the  OATS  model).  Also,  the 
target  is  assumed  to  be  stationary  (i.e.  the  change  in 


sin  i'P)  =  -  sin  (C) 


sin’^l^psin  (C)  j 
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FIGURE  2  RANGING  GEOMETRY 


distance  from  the  target  to  the  satellite  is  due  only  to  the 
velocity  of  the  satellite).  The  satellite  is  set  in  a 
circular  orbit  of  some  selectable  altitude  and  some  select¬ 
able  offset.  The  offset  is  simply  the  distance  that  the 
target  is  from  the  satellite’s  ground  track  at  its  closest 
point  of  approach.  The  offset  is  along  the  Y  axis  and  the 
satellite  ground  track  is  parallel  with  the  X  axis.  The  x-r 
plane  intersects  the  earth’s  surface  at  the  origin  only  and 
therefore  the  Z  axis  represents  a  normal  to  the  earth. 
Expanding  on  Figure  2,  the  following  terms  apply: 

d  -  distance  from  the  target  to  the  satellite 

y  -  shortest  distance  from  the  target  to  the  satellite's 
ground  track 

h  -  altitude  of  the  satellite  above  the  earth’s  surface 
d„  -  the  "no-offset  range";  defined  as  the  distance  from 

n 

the  target  to  the  satellite’s  projection  on  the 
plane;  this  is  the  shortest  distance  from  the 
satellite  to  the  Y  axis. 

-  zenith  angle  of  the  satellite  measured  from  the  2  a.'' i  s 

6  -  the  "no-offset  zenith  angle'  defined  as  the  angle 

from  Z-axis  to  d  vector 

The  known  data  entries  into  the  OATS  model  are  satellite 
altitude,  h,  the  offset,  y,  and  the  no-offset  zenith  angle. 

0  .  To  calculate  d,  the  following  equation  can  be  used 

n 

(9:385  )  : 


To  solve  for  d,  the  terms  x,  y  and  z  must  be  defined.  The 
satellite,  if  in  a  perfect  circular  orbit  over  a  perfect 
spherical  earth,  has  a  constant  altitude,  h.  In  Equation 
15,  2  is  the  apparent  height  of  the  satellite  above  the  x-y 
plane.  In  the  OATS  model,  z  must  always  be  positive  since 
the  no-offset  zenith  angle,  is  limited  from  80'’to  -60^. 
The  value  for  z  reaches  its  maximum  when  is  zero  and  it 
reaches  its  minimum  when  is  ±80  degrees.  The  value  for  z 
in  Equation  15  is  a  function  of  S  such  that 

N 

2  =  d„  cos(0^ )  (  1  C  I 

rl  n 

where  d^  is  the  no-offset  range.  The  apparent  satellite 
altitude  is  now  a  function  of  the  no-offset  range. 

The  value  for  y  is  not  straightforward.  Since  tne 
earth  is  curved,  y  does  not  exactly  equal  the  ground  range 
of  the  offset.  If  y  =  d^  is  the  "direct  orthogonal  range  , 
then  let  D  be  the  ground  range.  Figure  3  describes  the 
geometry.  From  the  figure. 


and  the  circumference  of  the  earth  is  2i!R  .  The  portion  of 
the  circumference  that  forms  the  angle  is 


direct 

orthogonal 

distance 


or  <PR 


t  18; 


In  the  figure  therefore,  D  =  Still  assuming  that  the 

zenith  angle  is  ±80^’  then  the  distance  d  is  strictly  a 

O 

function  of  D  such  that 


d  =  R  (ain't*) 

o  e 


(  1  9  I 


For  small  offset  distances,  d^  «  0.  The  error  in  this 
approximation  is  given  in  percentage  by 


%e 


do-D 

t) 


X  100 


(  20  ) 


Figure  4  graphically  portrays  this  error.  Since  the  error 
is  so  small  out  to  500  km,  the  model  uses  the  approximation 
d^  =  D,  the  direct  orthogonal  range  is  approximated  by  the 
ground  range.  For  Equation  15  then,  y  =  d^  and  the  OATS 
model  is  limited  to  offsets  of  500  km  or  less.  Beyond  500 
km,  the  user  must  accept  the  offset  error  or  compensate  with 
some  other  means. 

The  final  term  in  Equation  15  that  must  be  identified 
is  X.  Using  the  same  logic  that  formed  the  basis  for 
Equation  16,  the  value  for  x,  the  distance  out  the  X  axis  is 


X  =  d^  cos(0^) 


(21  ) 
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Equation  21  shows  that  x  Is  a  function  of  the  zero  offset 
range,  d„.  The  planar  geometry  for  this  case  (i.e.  d  =  d^ 
since  y  =  0)  is  shown  in  Figure  5.  The  range  from  the 
target  to  the  satellite  can  be  approximated  with  Equation  22 
shown  below. 


A  A  h 

Therefore, 

**  cos  ©ij 


(22) 


(23) 


This  approximation  is  only  good  for  small  angles  6^  because 
it  assumes  no  curvature  in  the  satellite’s  path. 

To  be  more  precise  and  take  curvature  into  account,  the 
law  of  cosines  can  be  used  to  precisely  calculate  d^^.  This 
algorithm  for  the  computation  of  the  distance  d^,  was 
developed  with  Don  Holland  and  a  version  of  it  is  listed 
below  (11).  A  general  form  of  the  law  of  cosines  is  listed 
in  Equation  24  (9,0). 

c^  =  a^  b^  -  2  a  b  cos(C)  (24) 


Applying  Equation  24  to  Figure  5  yields 

r^  =  +  df  -  2  R  d^cos(n  -  9^)  (25) 

•  If  e  If  N 
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SATELLITE 


FI6URE  5  NO  OFFSET  PLANAR  GEOMETRY 
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Since  r  =  R  +  h  and  R  ,  h  and  0  are  entering  arguments 

e  e  II 

Equation  25  can  be  rewritten  as 


d^  -  2  R  d„cos(jr  -  0^)  +  R^  -  (R  +  h)^  =  0 

N  e  N  Nee 


letting 


then 


K,  =  -2  R  cos(Tr  -  0  ) 

3  e  N 

and 

K  =  R^  -  (R  +  h)^ 

^  e  e 


^  *<3  ^  =  0 


(  26  ) 
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(  29  ) 


Solving  for  d„  with  the  quadratic  equation 

N 


_  -K,  ±  /  -  4 


(  30  ) 


Since  0^  can  range  from  o'’  to  90°  (lower  than  the  horizontal 
is  not  allowed  in  the  OATS  model),  the  term  will  always 
be  positive.  For  all  satellites  with  a  positive  altitude, 
the  term  will  always  be  negative.  Therefore,  the  term 
under  the  radical  will  always  be  a  positive  quantity. 
Equation  30  can  be  rewritten  as 


-K  +  /  kI-  i  K 

H  - T - 


(31 
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Therefore,  Equation  31  always  returns  a  real  and  positive 
value  for  the  range  from  the  target  to  the  satellite. 

Recall  d  is  the  zero  offset  range  from  satellite  to 
target. 

Going  back  to  Figure  2,  d^^  is  now  simply  a  function  of 
satellite  altitude,  h  ,  radius  of  the  earth,  R  ,  and  the 

e 

no-offset  zenith  angle,  The  term  becomes  a  key 

component  in  the  OATS  model. 

The  three  key  terms  of  Equation  15  (x.y,z)  have  been 
defined  with  Equations  16,  21  and  the  y  =  d  approximation. 

o 

The  distance  from  the  satellite  to  the  target  can  therefore 
be  computed  with  good  accuracy  out  to  1000  km  offset  range 
and  beyond  1000  km  offset  if  Equation  19  is  used  to  correct 
long  range  variation  in  y. 

It  must  be  noted  however  that  this  range,  d,  is  the 
geometric  range  and  is  not  the  optical  path  range.  Tne 
actual  optical  path  through  the  atmosphere  would  be  greater 
than  or  equal  to  the  direct  range,  d,  due  to  the  refraction 
of  the  light  through  the  atmosphere. 


The  process  of  computing  the  transmittance  of  the  laser 


energy  through  the  atmosphere  was  very  labor  intensive  and 
yet  yielded  the  model’s  least  flexible  results.  The  labor 
intensive  portion  of  the  transmittance  calculations  was  from 
using  FASC0D2  on  the  Aeronautical  Systems  Division’s  CDC_ 
Cyber  computer  system.  The  lack  of  flexibility  in  the 
results  stems  from  the  narrow  applicability  of  the  results. 
Over  120  FASCODE  computer  runs  yielded  only  two  types  of 
very  specifically  defined  atmosphere  types. 

Transmittance  in  this  model  is  simply  the  ratio  of  the 
signal  that  gets  through  the  atmosphere  to  the  signal  that 
was  sent  from  the  radiating  source.  The  transmittance,  i. 
therefore  varies  from  0  (no  power  transferred)  to  1  (all 
power  transferred).  As  previously  discussed,  transmittance 
through  the  atmosphere  is  a  function  of  wavelength.  This 
wavelength  dependence  is  determined  by  the  composition  of 
the  atmosphere. 

FASC0D2  is  constructed  so  as  to  give  the  modeler  a  wide 
combination  of  atmospheres  to  simulate.  Thousands  of 
combinations  are  available.  These  combinations  are  due  to 
the  ability  of  the  FASC0D2  user  to  define  a  wide  number  of 
constituents  in  the  atmosphere  at  different  altitudes  and 
concentrations.  Items  that  can  be  selected  in  various 
concentrations  include,  up  to  28  molecular  species,  clouds, 
rain  and  even  volcanic  particles.  Six  pre-established 
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atmospheric  models  are  available  to  define  general  atmos- 

I  phere  profiles  for  certain  conditions  (1:1).  One  of  these 

pre-established  models,  the  Midlatitude  Winter  model  was 
used  to  generate  the  transmittance  data  for  the  OATS  model. 

I  This  model  is  established  for  latitudes  about  45°  North 

Latitude  (1:1).  It  therefore  covers  much  of  the  United 
States,  Europe  and  Asia,  where  aircraft  operations  are 
common.  The  Mid  latitude  Winter  model  defines  the  temper¬ 
ature  profile  of  the  atmosphere  which  is  used  with  a  cor¬ 
responding  data  base  defining  the  concentrations  of  gases  at 
different  altitudes  (1,1-3). 

The  variable  IHA2E  was  the  only  variable  changed  in 
order  to  create  two  types  of  atmospheres  for  the  OATS  model . 

1  The  first  atmosphere  had  IHAZE  held  constant  at  0,  defining 

no  haze  in  the  atmosphere.  The  second  atmosphere  had  IHAZE 
held  to  2,  defining  a  5  km  rural  haze  profile.  Although 

I 

further  runs  were  attempted  using  clouds,  unresolved 
software  anomalies  on  the  CDC  cyber  computer  prevented 
success  with  cloud  models. 

Each  run  of  FASCODE  was  designated  at  a  set  zenith 
angle  and  resulted  in  a  series  of  17  transmittance  values 
corresponding  to  17  wavenumbers  or  wavelengths.  The  spacing 
between  each  increment  of  wavenumber  was  approximately  x 
where  x  =  (2.4  x  10'^  ♦  wavenumber),  an  extremely  narrow 
section  of  the  electromagnetic  spectrum  at  laser  frequen¬ 
cies.  By  choosing  the  center  frequency  transmittance  values 
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for  runs  at  a  variety  of  zenith  angles,  a  curve  was  estab¬ 
lished.  Eight  curves  were  developed,  one  for  each 
atmosphere  and  laser  combination.  The  curves  defined  the 
relationship  of  transmittance  to  zenith  angle. 

In  order  to  be  able  to  generate  transmittance  values 
for  any  zenith  angle,  an  equation  for  each  curve  was  sought. 
SAS^  (21)  software  was  used  for  this  purpose.  Using  the 
General  Linear  Model  (GLM)  procedure,  coefficients  for  a 
fifth-order  polynomial  equation  were  computed  for  each 
curve.  This  procedure  uses  a  least  squares  fitting 
algorithm  in  its  computations  (21:62). 

The  precision  of  the  curve  fit  was  enhanced  through  a 
data  transformation.  Instead  of  using  the  zenith  angle  as 
the  independent  variable,  the  cosine  of  the  zenith  angle  was 
used.  The  input  data  for  this  process  is  listed  in  Appendix 
C  for  each  curve. 

The  accuracy  of  the  curve  fit  with  the  actual  data  was 

very  good.  A  means  to  determine  the  goodness  of  fit  is  the 
2 

R  value  that  is  automatically  generated  for  each  GLM  run. 

2 

"R  measures  how  much  variation  in  the  dependent  variable 
can  be  accounted  for  by  the  model."  (21:64)  An  R^  value 
better  than  ,999  was  achieved  for  each  of  the  eight  curves. 
In  other  words,  the  fifth-order  polynomial  equation  was 
almost  a  perfect  fit  to  the  curve  in  all  eight  cases. 

Data  input  for  the  curve  fitting  were  restricted  to 
zenith  angle  ranges  from  o'*  to  a  maximum  of  80°.  Although 
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some  data  was  generated  with  FASCODE  beyond  the  limit  of 
80^,  this  data  was  not  used.  For  example,  the  CO^  with 
carbon-12  curve  with  no  haze  was  developed  for  zenith  angles 
from  o“  to  90°.  This  particular  series  showed  a  transmit¬ 
tance  curve  that  appeared  exponential  for  the  first  80°  but 
had  a  tail  in  the  last  10°  (i.e.  from  80°  to  90°).  In 
accounting  for  this  tail,  the  curve  fitting  algorithm  drove 
the  accuracy  of  the  first  80°  off  by  a  small  margin.  The 
data  input  was  subsequently  restricted  to  a  maximum  of  80° 
and  resulted  in  an  improved  accuracy  in  the  fit.  This 
increase  in  accuracy  was  traded  for  the  extended  range  (  30° 
to  90°)  for  the  OATS  model. 

The  80°  upper  limit  also  makes  sense  in  a  realist''c 
scenario  of  a  laser  radar  on  a  satellite  tracking  low 
altitude  targets.  Further  complications  to  the  system  arise 
in  80°  to  90°  sector.  The  additional  attenuation  in  this 
regime  (not  modeled  with  FASCODE  here)  would  be  very  high 
from  dust,  pollen,  terrain,  etc.  for  low  altitude  targets. 
Also,  refraction  would  be  at  its  worst. 

With  the  coefficients  from  the  polynomial  curve  fit, 
one  can  easily  calculate  a  very  accurate  transmittance  for 
any  zenith  angle.  However,  the  underlying  assumption  here 
is  that  FASCODE  is  a  perfect  model  of  the  atmospheric 
transmittance  and  that  the  atmosphere  exactly  matches  the 
atmosphere  loaded  into  FASCODE. 

Because  the  atmospheres  modeled  by  FASCODE  (  i.e., clear 
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atmosphere  with  no  aerosols  or  5  km  haze  the  only  aerosol) 
yield  relatively  optimistic  values  of  transmittance,  an 
additional  attenuation  capability  was  built  into  the  OATS 
model.  Although  unsubstantiated  for  accuracy,  an  exponen¬ 
tial  model  was  used  to  simulate  the  variance  of  the  trans¬ 
mittance,  T,  as  a  function  of  zenith  angle.  The  value 
entered  into  the  model  is  t.he  additional  attenuation  desired 
for  a  zenith  angle  of  zero  (in  dB).  For  angles  beyond  0 
degrees,  the  additional  attenuation  is  computed  with  the 
following  equation 


where  L  is  the  length  through  the  atmosphere.  For  this 
equation,  the  atmosphere  is  assumed  to  end  at  an  altitude  of 
120  km,  the  definition  of  space  used  for  the  FASCODE  runs. 
The  value  for  L  is  therefore  at  a  minimum  when  0  is  zero 
(L=120  km)  and  at  a  maximum  when  9  is  at  80°,  the  OATS  limit 
(L'»557  km).  The  geometry  is  therefore  the  same  as  that 
described  in  the  ranging  algorithm.  Since  L  is  a  function 
of  6  only,  the  value  for  L  is  simply  the  range  if  h  is  120 
km  for  all  values  of  9.  Therefore,  the  ranging  function  is 
invoked  in  the  OATS  model  to  calculate  the  value  L  for  any 
zenith  angle. 

The  coefficient  a  is  determined  from  the  input 
attenuation  requirement,  5^(0),  The  term  3  is  the 
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attenuation  input  by  the  modeler  and  3(0)  designates  the 
attenuation  at  a  zenith  angle  of  zero  degrees.  Since  3(6) 
represents  the  attenuation  in  dB,  the  ratio  of  received 
power  to  transmitted  power  is. 


lower 


Tx  power 


where  n  ~ 


-3(0) 

10 


(33) 


so  that  the  attenuation  coefficient  a  can  be  defined  as 


(34) 


Therefore,  the  transmittance  factor,  is  defined  as 


xJ6)  =  e 


-aL(e) 


(35) 


The  transmittance  factor,  is  multiplied  by  the 
computed  transmittance  from  the  fifth-order  polynomial  for  a 
clear  atmosphere  to  yield  the  overall  model  transmittance  at 
a  given  zenith  angle. 

The  transmittance  through  the  atmosphere  can  be 
calculated  in  the  OATS  model  in  three  ways: 

1)  a  clear,  no-aerosol  atmosphere  made  up  of  28  molecular 
species 

2)  atmosphere  number  one  with  5  km  rural  haze 

3)  atmosphere  number  one  with  a  selectable  additional 
attenuation  factor. 
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Lidar  Cross  Section 

One  of  the  more  complicated  portions  of  the  OATS  model 
is  the  generation  of  Lidar  Cross  Section  (LCS)  values  for 
the  targets  that  the  lidar  system  is  intended  to  track.  The 
LCS  of  a  target,  *t,  is  analogous  to  the  Radar  Cross  Section 
(RCS)  of  a  target  in  the  microwave  frequency  range  but  is 
merely  differentiated  in  its  application  to  the  optical 
frequency  ranges. 

Bachman,  in  his  discussion  of  optical  cross  sections  of 
targets,  starts  his  analysis  with  the  equation 

o  =  p  G  A  (36) 


'■ .  .  .where 

p  =  reflectance  of  the  target  surface 
G  =  gain  of  target 

A  =  projected  physical  area."  (2:35) 
He  goes  on  further  to  define  gain,  G  as 


(37) 


"where  A^  is  defined  as  area  of  coherence  (equivalent  normal 
incidence  flat  plate  area)  of  the  target  and  0^  is  the 
backscatter  solid  beamwidth  from  the  target."  (2:35) 
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Bachman  continues  his  development  and  says,  "...the 
equivalent  specular  scatterer  is  equivalent  to 


which  is  identical  to  the  radar  cross  section  (RCS)  at 
wavelength  X."  (2:36)  Indeed,  Skolnik  cites  Mentzer  and 
defines  the  RCS  of  a  "large  flat  plate  of  arbitrary  shape" 
(23:43)  with  the  same  equation. 

However,  Bachmart  further  goes  on  to  define  the  cross 
section  of  a  diffuse  scatterer  by  defining  =  rt  and 
stating  that  the  relationship  "assumes  backscatter  energy 
over  n  steradians  of  the  4n  steradians  of  the  isotropic 
sphere..."  (2:36)  With  this  assumption,  he  ends  up  with  the 
cross  section  of  a  diffuse  reflector  defined  as 

a  =  4  p  A  COS^6  ( 39  ) 

In  the  OATS  model,  two  types  of  targets  are  available, 

the  flying  sphere  and  the  flying  brick.  The  spherical 

target  approximation  is  the  simplest  in  its  derivations  and 

application.  From  Skolnik,  who  cites  Mentzer,  the  cross 

2 

section  of  a  sphere  is  na  where  a  is  the  radius  of  the 
sphere  (23:43).  However,  in  this  treatment  of  cross 
section,  the  target  is  assumed  to  be  Lambertian  or  a  diffuse 
reflector.  Therefore,  based  on  Equation  39,  one  can  show 
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that  the  ideal  diffuse  reflecting  sphere  has  a  LCS  of 


-|-  n  p  (40) 

where  r  is  the  radius  of  the  sphere. 

This  can  be  done  by  integrating  over  a  hemisphere  so 

that 

P»r/2 

a  =  4  p  A  cos^Q  dO  (41  ) 

OiO 

Following  the  integration  scheme  used  by  Slater  to  solve  a 
similar  problem,  the  increment  of  integration  is  a  ring 
representing  a  segment  of  the  hemisphere  defined  by  d©. 
According  to  Slater’s  orientation  in  this  type  of 
integration,  the  area  of  the  ring  (if  radius  =  1)  is 

A=2nsin(6)d6  (42) 

(24:529) . 

Combining  Equations  41  and  42  yields 

71/2 

o  =  8  n  p  r  ^  sin(6)  cos^(0)  60  (43) 

eio 

The  solution  to  this  equation  is 

o  -  -|-  n  p  (44) 
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The  flying  brick  model  uses  three  orthogonal  surfaces 
oriented  in  a  manner  such  that  no  surface  overlaps  or 
shadows  the  other  (e.g.  a  brick).  This  model  was  developed 
to  more  accurately  represent  an  aircraft.  Although  most 
aircraft  are  not  shaped  like  a  brick,  most  aircraft  tend  to 
have  a  larger  projection  area  (silhouette)  from  a  top 
perspective  than  from  a  frontal  perspective.  By  making  the 
assumption  that  a  larger  projection  area  equates  to  a  larger 
LCS,  one  can  assume  that  most  aircraft  have  a  larger  LCS 
from  a  top  perspective  than  from  a  frontal /rear  perspective. 
Since  a  satellite-borne  lidar  system  would  potentially  track 
a  target  aircraft  over  a  160°  arc  (the  OATS  model  limit) 
through  the  aircraft’s  upper  hemisphere,  a  varying  LCS  would 
be  a  reasonable  assumption.  The  original  idea  for  a  brick 
came  from  a  discussion  with  Lt.  Blake  and  Mr.  Weigand  (27). 
The  brick  model  provides  the  changing  LCS  that  could  be 
expected  from  a  true  aircraft  as  the  satellite  transits 
through  a  full  range  of  zenith  angles. 

If  the  aircraft’s  LCS  is  known  and  can  be  simplified  to 
a  sum  of  three  perpendicular  cross  sections,  those  data 
could  be  put  into  the  OATS  model  inputs  as  the  area  and 
reflectivity  of  the  target. 

The  following  assumptions  form  a  foundation  for  the 
OATS  model : 

1.  The  frontal  and  rear  faces  of  the  brick  have  equal 

projected  areas. 
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2.  The  side  faces  of  the  brick  have  equal  projected  areas. 


3.  The  aircraft  (the  brick)  is  in  straight  and  level  flight 
with  no  crab  (i.e.  the  true  heading  vector  is  normal 
to  the  frontal  face). 


4.  The  skin  reflectance,  p,  is  uniform  and  the  surface  is 
diffuse  (or  Lambertian). 


5.  The  sum  of  the  cross  sections  of  the  three  surfaces  of 
the  brick  is  equal  to  the  total  cross  section. 

This  geometry  can  be  pictured  in  Figure  6.  The  frontal  or 
rear  projected  areas  are  depicted  with  PA1 ,  the  top  pro¬ 
jected  area  as  PA2  and  the  side  projected  area  as  PA3.  The 
satellite  ground  track  and  target  geometry  are  further 
depicted  in  Figure  7.  To  maintain  a  consistent  three  face 
object,  a  geometry  transformation  was  developed.  This 
transformation  realigns  the  target  faces  so  that  the  normal 
vector  to  the  frontal  face  is  always  parallel  with  the 
satellite’s  geometry  (see  Figure  8  and  9).  From  the  new 
configuration,  symmetry  about  the  2  axis  can  be  assumed  so 
that  the  variation  in  offset  zenith  angle,  from  90°  to 
0°  is  symmetrical  with  the  variation  from  0°  to  -90°. 

If  a  simple  projection  area  transformation  is  used,  the 
effect  is  obscured.  The  simple  transformation  is 


s  = 

IN 

cos 

< 

+ 

sin  P 

(45) 

s  = 

In 

cos  ^| 

< 

+ 

sin  p 

(46) 

These  transformations  simply  change  the  projected  areas  of 
the  target.  The  top  projected  area  remains  the  same, 


iurfacec  are  orthogonal  with  each  other) 


FIGURE 


FIGURE  7  ACTUAL  TARGET  ORIENTATION 


TAR6ET 


FIGURE  9  TRANSFORMED  TARGET  GEOMETRY 


independent  of  the  targets  true  heading  if  the  target  is 
straight  and  level.  However,  the  failure  of  these 
transformations  is  found  in  the  fact  that  the  variation  in 
planar  LCS  is  a  function  of  cos^d  (from  Equation  39). 
Therefore  the  sum  of  projected  areas  in  Equations  45  and  46, 
when  viewed  from  an  angle  does  not  have  the  same  LCS  as 
the  transformed  areas  (A’^and  A’^)  due  to  the  cos^B  term. 

The  following  example  illustrates  this  point  and  the 
transformation  inaccuracies.  The  true  frontal  projected 
area  is  5  square  meters,  the  side  projected  area  is  20 
square  meters.  The  angle  B  is  30”.  Therefore, 


a;  = 

1 5  cos  ^1  +  20 

sin  ^ 

«  14.33 

sq 

meters 

(47  ) 

^3  = 

|20  cos  /9(  +  5 

sin  /? 

19.82 

sq 

meters 

(48  ) 

To  simplify  the  example,  the  zenith  angle,  0,  is  set  to  90\ 
Note,  the  top  surface  does  not  contribute  if  =  90”.  The 
target’s  untransformed  cross  section  is  ‘^3  where 

<1,  =  4  p  (  5  )  cos^S  (  49  ) 

=  4  p  (20)  cos^S  (50) 

If  the  angle  5  is  set  to  45”  for  this  example  (where  <3 
defines  the  satellite’s  position  along  its  path),  then 


n 


and 


<7^  =  10  P 
<73  =  40  p 


so  that  <7^  =  50  p 


But  the  transformed  geometry  yields 


<7^  **  28.66 
«  39.64 


r: 


i 


<7^  =  68.30  square  meters. 

The  transformations  in  Equation  45  and  46  are  therefore 
shown  to  be  imprecise. 

To  correct  for  this,  the  transformation  equations  were 

2 

modified  to  reflect  the  cos  0  dependence  to  be 


t 

N  = 

cos^Pl  +  A3  sin^^ 

(51  ) 

^3  = 

cos^/Sj  +  sin^P 

(  52  ) 

These  transformation  equations  will  accurately  transform  the 
surface  areas  to  pseudo-areas  such  that  the  psuedo-areas  can 
be  used  in  the  LCS  equations  for  LCS  computations.  The 
psuedo-areas  do  not  represent  true  projection  areas  because 
of  the  square  terms. 

It  is  important  to  point  out  that  although  the  accuracy 
of  the  transformation  is  improved,  the  entire  procedure 
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rests  on  the  premise  that  the  flying  object  consists  of 
three  orthogonal  surfaces.  As  the  three  surface  approxi¬ 
mation  varies  from  reality,  so  does  the  accuracy  of  the 
model.  The  three  surface  model  does,  however,  give  the 
modeler  the  opportunity  to  have  the  target  LCS  realistically 
change  with  aspect  angle. 

The  next  step  in  the  construction  of  the  flying  brick 
model  was  the  development  of  equations  to  accurately  compute 
the  individual  surface  LCS  values.  From  Equation  39,  it  is 
apparent  that  the  off-normal  aspect  angle  needs  to  be  cal¬ 
culated  or  estimated. 

The  simplest  case  involves  no  offset  (see  Figure  10), 

In  this  case,  only  two  surfaces  are  involved,  PA1  and  PA2, 
when  p  =  0  or  Tt.  Even  if  the  targets  heading  is  not  paral¬ 
lel  with  the  satellite’s  ground  track  (i.e.,  0^0  ,  0  ^  ri  )  , 
the  transformation  in  Equations  51  and  52  yield  an  equiva¬ 
lent  geometry  such  that  the  pseudo-areas  Aj  and  A^  can  be 
used  to  keep  the  geometry  simple,  i.e.,  only  two  side  sub¬ 
ject  to  illumination/reflection.  With  no  offset, 


o  =  a  +  0 

T  t  2 
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where 


=  4  p  A  cos^  ~ 

=  4  p  A  cos^6  (  55  ) 
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SATELLITE 


FIGURE  10  ZERO  OFFSET  GEOMETRY 


HOWEVER,  by  including  an  offset  into  the  geometry,  the 


model  switches  to  the  geometry  depicted  in  Figure  11. 

Using  the  angle  relationship  for  spherical  coordinates 
(9:385)  and  transferring  it  into  the  present  coordinate 
system  and  symbology. 


(56  ) 


where  angle  from  the  2  axis  to  the  vector  formed 

by  the  distance  from  the  target  to  the  satellite.  This  is 
the  off-normal  angle  for  PA2  (or  A^)  of  the  brick  model.  In 
terms  of  the  no-offset  zenith  angle, 


so  that 


cos  B 


Z  =  d,  cos  e. 


(  57  ) 


(58) 


where  d^  is  the  no-offset  range.  The  value  for  x  is 

n 


X  =  d^  sin  0^ 


(  59  ) 


Since  d„  is  already  computed  during  ranging,  Equation  58  and 

fl 

59  are  convenient.  Using  the  same  relationship  as  Equation 
56,  a  and  «  are  given  as  follows. 
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FIGURE  1 1  OFFSET  GEOMETRY 


X 


=  COS 


2^2  2 

X  +  y  +  z 


(60) 


=  cos 


2^2^  2 

X  +  y  +  z 


(61  ) 


From  these  relationships,  the  individual  surface  LCS  values 
can  be  computed.  The  following  equations  form  the  individ¬ 
ual  LCS  values  that  are  necessary. 


=  4  p  A’  cos 
=  4  p  A,  cos  a 


4  p  A^  COS 


The  total  brick  LCS  from  any  aspect  is  therefore 


■j 


<T  =  a  +  o  +  a 

X  1  2  3 


(  62  ) 


(  63  ) 


(  64  ) 
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The  Lidar  Range  Equation 


The  fifth  major  subsection  of  the  OATS  model  is  the 
lidar  range  equation,  its  form,  its  variables  and  the 
underlying  assumptions.  This  equation  is  the  hub  of  the 
OATS  model . 

The  particular  form  of  the  lidar  range  equation  used  in 
the  OATS  model  is  taken  from  Dr.  Paul  McManamon’s  analysis 
of  lidar  ranging  along  with  his  equations  for  transmission 
beam  gain,  G  ,  signal  to  noise  ratio,  S/N,  and  effective 
aperture,  A  (17:1).  These  equations  are  listed  below. 


P„  =  P^  ^ 

Tx  Tx  4  rt  ^  ^  p2 


A  = 

e 


S/N  = 


Prx  ^ 


( t>6  ) 


\  67  ) 


(  68  I 


The  prior  three  subsections  defined  the  target’s  lidar 
cross  section,  c»^,  the  atmospheric  transmittance,  t,  and  the 
range  from  the  target  to  the  satellite,  R,  for  use  with 
Equation  1.  The  following  pages  will  define  the  terms  from 
Equation  1  that  have  not  been  defined  or  explained,  specif¬ 
ically,  A^,  G^^,  P^^  and  P^^. 
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Equation  66  defines  the  receiver  aperture,  A.  This 

e 

equation  is  valid  for  a  circular  receiving  aperture  of 
diameter,  D,  and  the  aperture  is  simply  the  area  defined  by 
a  circular  optic  with  a  diameter,  D.  However,  Stimson,  in 
his  analysis  of  antenna  gain,  also  uses  and  says, 
"Effective  area  is  equal  to  physical  area  times  aperture 
efficiency"  (26:137)  where  he  later  labels  aperture 
efficiency,  h,  (hereafter  designated  n^).  Combining  both 
concepts,  the  new  equation  would  be 

n  0“  n 

A,  =  - 4 - 

The  aperture  efficiency  here  could  be  less  than  one  for 
reasons  such  as:  1)  an  occlusion  from  satellite  structure, 
2)  occlusion  from  optics  (like  a  secondary  cassegrain  re¬ 
flector)  and  3)  less  than  10056  of  the  circular  optic  being 
used  for  collection  (e.g.  structure  attachment  points  on  a 
parabolic  dish).  However,  the  OATS  model  uses  Equation  66 
as  it  is  written  thereby  assuming  an  aperture  efficiency  of 

one.  The  modeler  must  be  cognizant  of  this  and  if  n.  is 

* 

less  than  one,  the  value  for  0  must  be  altered  so  to  result 
in  a  properly  proportioned  value  for  A  . 

Equation  67  defines  transmitter  gain  from  beam  direc¬ 
tivity.  The  OATS  model  deviates  from  Equation  67  on  this 
term  and  uses  a  more  precise  circular  aperture  assumotion. 
The  new  relationship  is  listed  below  as  Equation  70. 


(  70  ) 
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In  Equations  67  and  70,  6  is  the  "transmit  beam  angular 
diameter"  (17:1)  or  the  full  angle  beam  divergence  measured 
at  the  3  dB  points  (17:2).  The  deviation  between  Equation 
67  and  Equation  70  results  from  the  basic  definition  of 
gain.  Stimson  defines  the  gain  due  to  the  directivity  of  an 
antenna  in  the  following  way, 

"The  gam  of  an  antenna  is  the  ratio  of  the  power  per 
unit  of  solid  angle  radiated  in  a  specific  direction 
to  the  power  per  unit  solid  angle  that  would  nave 
been  radiated  had  the  same  total  power  been  radiated 
uniformly  in  all  direction..."  (26:137) 

Or  rewriting  Skolnik’s  gain  definition  in  present  terms. 
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4  n 
Q 


i  7  1  ) 


where  the  term  0  is  defined  as  the  beam  area  by  Skolnik 
(23,262).  However,  for  a  circular  antenna,  0  is  a  function 
of  6  ,  the  beam  divergence  angle.  Using  Slater’s  definition 
of  solid  angle,  here  defined  as  ft, 

A 

0  =  -j-  (72) 

r 

S 


where  A  is  the  projected  area  on  a  sphere  of  radius  r 

8  C 

(24:528).  If  the  radius  of  the  sphere  (distance  from  the 
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beam  origin  to  projected  area  plane)  is  1,  then  0  =  A  , 
Assuming  r  =1,  for  small  divergence  angles, 

8  D 

so  that  6^  ^  d  ( 73  ) 


Y-  tan 


where  ct  is  the  diameter  of  the  projected  circular  area. 

2 

Since  the  projected  area  is  n  r^  where  d  =  2  r^  ,  then 


Q  = 

SA  4  B 


I  74  ) 


SO  that  the  gain  of  a  circular  antenna  is 


4  n 


16 


'tx  "  Q. 


(  75  ) 
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The  OATS  model  uses  this  relationship  and  is  therefore  built 
on  the  assumption  of  a  circular  beam. 

Notice  that  the  gain  of  the  antenna,  is  based  on 

the  beam  divergence,  of  the  laser  where  is  definea  bv 
the  3dB  points.  This  is  an  optimistic  approximation  since 
only  one-half  of  the  laser  power  is  within  the  angle  defined 
by 

The  final  terms  requiring  explanation  in  Equation  1  are 
the  two  power  terms,  and  .  These  terms  are  simply 
the  power  transmitted,  P^  ,  and  the  power  received,  P„  . 

Tx  Px 


The  transmitted  power  is  the  power  in  a  single  laser  pulse 


and  can  be  calculated  by  dividing  the  energy  per  pulse  by 
the  duration  of  the  pulse  (pulse  width).  The  received  power 
is  the  amount  of  energy  per  time  period  that  has  reflected 
from  the  target  and  is  collected  in  the  receiver’s  optics. 

To  add  more  flexibility  into  the  model  a  processing 
gain  factor,  <5,  was  included.  This  factor  is  nonspecific  in 
the  OATS  model  and  can  be  used  to  apply  the  additional 
increase  in  the  signal  from  any  post  reception  processing 
technique.  It  allows  the  modeler  to  include  gains  in  signal 
strength  from  such  techniques  as  coherent  integration  of  tne 
returned  pulses  but  is  flexible  enough  that  it  can  be  usea 
for  any  beneficial  or  detrimental  designs  in  the  system. 

The  application  of  this  factor  is  therefore  very  design 
dependent. 

The  final  portion  of  the  lidar  range  equation  that  must 
be  addressed  is  the  signal-to-noise  ratio,  S/N.  The  OATS 
model  assumes  a  shot  noise-limited  sensor  and  a  bandwiath- 
pulsewidth  product  equal  to  one  (17:2).  (For  example, 
assuming  a  10  MHz  bandwidth  for  the  receiver,  the  matcmng 
pulsewidth  would  be  100  nanoseconds  for  a  BW-PW  product  of 
one.)  The  assumption  of  a  shot  noise- 1 imi ted  sensor  in  the 
OATS  model  greatly  simplifies  the  process  of  predicting  the 
noise  in  the  system.  It  is  also  validated  by  Boyd  from  ms 
statement ,"... the  predominant  noise  source  in  a  well- 
designed  optical  heterodyne  detection  system  is  shot  noise 
in  the  photocurrent.  "  (4,204) 
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From  Equation  68,  B  is  the  bandwidth  (17,2)  or  more 

specifically,  the  "i-f  bandwidth"  (4,204)  which  is  the 

intermediate  frequency  bandwidth.  The  term  n  is  defined  as 

the  detector  quantum  efficiency  (4,128),  not  to  be  confused 

with  the  aperture  efficiency. 

The  quantum  efficiency,  h,  of  a  photoemi ssi ve  surface 
measures  its  performance  and  is  defined  as  the  number 
of  electrons  emitted  per  incident  photon^  The  quant^um 
efficiency  varies  approximately  from  10'  to  6  x  10' 
for  most  photocathodes  and  it  varies  as  a  function  of 
the  incident  radiation  wavelength.  (22:31) 

The  OATS  model  sets  the  quantum  efficiency  to  0.1, 
eliminating  the  details  of  individual  detectors.  However, 
this  is  another  function  that  could  be  combined  in  the 
processing  gain,  if  the  modeler  desired  to  vary  the 
quantum  efficiency  of  the  system.  For  example,  if  the 
actual  quantum  efficiency  to  be  modeled  was  0.2,  the 
processing  gain,  could  be  increased  by  2  dB  to  allow  for 
this  better  quantum  efficiency. 

From  Equation  68,  the  term  i»,  is  the  frequency,  ^  is 
Planck’s  constant  and  P„  is  the  power  at  the  receiver 
(defined  with  Equation  1).  To  achieve  a  BW-PW  product  of 
one,  the  bandwidth  is  pre-established  in  the  model  once  the 
pulse  energy,  E^^,  and  pulsewidth,  PW,  are  set.  One  can 
combine  Equations  1  and  68  (17:2)  which  yields  the  following 
set  of  relationships. 
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Simplifying  with  PW*B  =  1  and  h  =  0.1  then. 


S/N  = 


E 

Tx  t 


10  A  IV 


i  /  a  ) 


which  is  then  simplified  to  the  core  lidar  range  eauation 
for  the  OATS  model : 


S/N  = 


E_  <5  <r  X 

T  X _ I _ 

(  2.4979x10'‘^)  Og  r"* 
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Another  feature  that  was  inserted  in  the  comoutational 
equations  in  the  OATS  model  was  for  spot  size.  With  a  low 
flying  satellite  and  a  large  aircraft,  one  can  conceive  of 
the  scenario  when  the  spot  size  of  the  laser  is  smaller  that 
the  size  of  the  aircraft.  If  the  projected  area  of  the 
target  is  larger  than  the  spotsize,  then  the  lidar  range 
equation  needs  to  be  adjusted.  To  accommodate  this,  the 
OATS  model  compares  the  computed  spot  size  with  the  target's 
projected  area.  If  the  projected  area  of  the  target  is 
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larger  than  the  spot  size,  then  the  spot  size  to  projected 
I  area  ratio  is  computed,  x.  This  ratio  is  directly  inserted 

into  the  range  equation  to  correct  for  small  spot  size 
errors  induced  from  the  LCS  equations.  This  direct  inser- 
,  tion  is  reasonable  and  simple.  The  reasonability  stems  from 

the  linear  effect  of  projected  area  on  the  LCS.  The  sim¬ 
plicity  results  from  the  direct  insertion  as  opposed  to 
recalculation  of  the  LCS  for  the  new  areas.  It  is  important 
to  note,  however,  that  the  laser  energy  is  assumed  to  be 
proportionately  spread  among  all  three  surfaces  of  the 
target  (e.g.,  if  the  top  surface  of  the  brick  represents 
7056  of  the  projected  area  of  the  brick  then  70%  of  the  laser 
irradiance  is  assumed  to  fall  on  the  top  surface).  The 
actual  distribution  of  the  irradiance  would  be  complex  and  a 
function  of  the  pointing  error  probability,  brick 
orientation  and  spot  size. 

This  analysis  and  implementation  of  the  spot  size 
correction  is  also  applicable  to  the  sphere.  With  a  sphere, 
however,  as  x  gets  small,  the  equation  for  LCS  of  a  sphere 
approaches  the  equation  for  LCS  of  a  planar  surface  with  no 
off  normal  angle  (assuming  perfect  pointing).  The  LCS  under 
perfect  pointing  conditions  could  be  readily  calculated  by 
integrating  over  a  partial  sphere.  But,  since  there  is  no 
such  thing  as  a  perfect  pointing  system  (especially  at  these 
ranges)  x  is  directly  inserted  into  the  lidar  range  equation 
just  as  it  is  for  the  brick  model. 
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The  computation  of  the  spotsize  uses  the  3  dB  beam 

B 

divergence  angle.  The  spotsize  is  therefore  the  3  dB  spot- 
size.  This  represents  only  an  approximation  to  the  actual 
irradiance  pattern  at  the  target.  A  specific  irradiance 
pattern  would  depend  on  the  precise  optical  configuration  of 
the  laser  and  the  perturbations  the  beam  encountered  in  the 
atmosphere. 

For  both  the  brick  and  sphere  models,  the  x  insertion 
technique  is  an  approximation.  Actual  observations  should 
validate  these  assumptions  since  the  beam  pointing  error  and 
beam  jitter  would  be  distributed  about  the  perfect  pointing 
and  no  jitter  points  (assuming  no  bias).  These  assump¬ 
tions  therefore  are  built  around  the  mean  of  naturally 
occurring  random  pointing  errors. 

This  spot  size  correction  reveals  another  significant 
limitation  and  assumption  of  the  OATS  model.  To  use  either 
computational  method  for  calculating  LCS  (i.e.,  spot  greater 
than  projected  area  or  spot  less  than  projected  area),  the 
underlying  assumption  is  that  the  beam  is  perfectly 
geometrically  coupled  (i.e.,  spot  pattern  matches  target 
pattern).  In  other  words,  the  beam  spot  on  the  target  is 
such  that  it  completely  engulfs  the  target  or  it  illuminates 
a  portion  of  the  target  with  no  portion  of  the  beam  missing 
the  target.  This  essentially  boils  down  to  the  idealistic 
assumption  that  there  is  no  spillage  of  the  beam  energy, 
even  with  perfect  pointing.  For  example,  if  the  spot 


60 


pattern  is  circular  with  an  area  equal  to  n  but  the  target 
pattern  is  that  of  an  aircraft  but  also  of  area  tr,  there  is 
no  way  that  lOOSi  of  the  energy  will  illuminate  the  target 
since  their  patterns  are  not  matched. 

If  the  pointing  is  perfect  and  if  the  beam  spot  size  is 
large  compared  to  the  size  of  the  target  (floodlight),  then 
the  geometrically  coupled  illumination  assumption  is  not 
bad.  If  the  area  of  the  spot  is  on  the  order  of  the 
projected  area  of  the  aircraft,  significant  portions  of  the 
beam  may  not  be  on  the  target  due  to  pattern  mismatching. 

If  the  spot  size  is  very  small  compared  to  the  dimensions  of 
the  target  aircraft  (i.e.,  extremely  small  B  )  then  the 

H 

pointing  accuracy  will  tend  be  the  dominant  source  of  error 
(either  on  target  or  off  target). 

Equation  79  is  therefore  rewritten  with  x  inserted  to 
form  Equation  80. 


S/N 


E,  <S  a  X  X 

T  X _ t _ 

(2.4979x10*^'"')  ef 

D 


(80) 


The  Inputs 

The  portion  of  the  OATS  model  in  which  the  modeler  is 
most  involved  is  the  data  input  segment.  It  is  in  this 
segment  that  the  modeler  describes  the  scenario  for 
anal ysi s . 
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Table  1  lists  all  the  cued  inputs  for  the  OATS  model, 
some  of  which  are  displayed  only  if  the  prior  data  entries 
require  the  additional  data  fields. 

Certain  fields  in  the  input  segment  require  further 
description.  The  BEAMOD  field  cues  the  user  for  the  beam 
model  to  be  used  for  the  analysis.  It  allows  three  options. 
The  first  option  asks  for  the  desired  beam  width  (BMDVG)  in 
radians.  With  the  second  option,  the  beam  divergence  is 
assumed  to  be  diffraction  limited  and  uses  Equation  81. 

This  equation  is  used  to  compute  the  "SdB  full  width  dif¬ 
fraction  limited  beam  divergence."  (17:2)  This  is  defined 
below. 


I  =  1 .08  {  81  ) 

Therefore,  the  3  dB  diffraction  limited  BEAMOD  entry  cues 
®  the  model  to  calculate  the  beam  divergence  angle  from  the 

selected  effective  receiver  diameter  (assumes  receiver 
optics  are  used  for  beam  transmission  also).  This  likewise 
assumes  that  the  aperture  efficiency,  is  1.0  for  the 
transmission  of  the  beam. 

If  BEAMOD  option  3  is  chosen,  the  OATS  model  computes  a 
beam  divergence  using  Equation  82  (taken  from  17:2)  for  a 
beam  that  is  "not  quite  diffraction  limited".  (17:2) 

r 

0g  =  1 .30  (82  ) 
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The  next  segment  that  requires  additional  data  entry  is 
atmosphere  selection  option,  ATMOS.  For  atmosphere  model 
three,  the  user  is  cued  to  enter  an  additional  amount  of 
attenuation  (  in  dB)  for  the  zero  degree  zenith  angle 
geometry.  At  all  other  zenith  angles,  the  OATS  model  com¬ 
putes  the  additional  attenuation  to  add  to  the  clear 
atmosphere  model  (discussed  previously  under  Transmittance). 

The  next  data  field  requiring  further  data  is 
LCSMOOEL.  If  the  spherical  model  is  chosen,  the  OATS  model 
cues  the  modeler  for  the  radius  of  the  sphere.  If  the 
flying  brick  model  is  chosen  the  model  cues  the  modeler  for 
the  area  of  the  three  sides  of  the  brick.  These  areas  are 
entered  in  square  meters.  The  top  area  should  represent  the 
projected  area  of  the  target  from  a  top-down  perspective. 

The  frontal  area  (also  the  rear  view  area)  is  the  area  from 
a  front-on  perspective.  The  side  area  is  likewise  the  area 
of  the  target  as  perceived  from  a  side  perspective. 

The  Output  Data  and  Calculations 

The  OATS  model  provides  a  variety  of  output  data  forms 
for  the  modeler  to  use  for  analysis.  The  method  of  output 
data  selection  is  through  the  DATASLCT  field.  Three  options 
are  currently  available  for  data  selection; 


■ 
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1)  the  no-offset  zenith  angle  where  the  S/N  level  first 
meets  the  previously  entered  required  S/N 

2)  the  no-offset  zenith  angle  where  the  minimum  detectable 
signal  is  first  met 

3)  the  signal  power  and  S/N  for  a  variation  in  one  of  the 
nine  selectable  variable. 

The  first  of  the  three  options  outputs  a  no-offset 
zenith  angle.  This  angle  represents  the  point  in  the 
satellite’s  orbit  at  which  the  satellite  can  distinguish  the 
signal  return  from  the  shot  noise.  Di sti ngui shabi 1 i ty  is 
established  by  the  SNdB  (required  S/N)  field  in  the  input. 
For  example,  if  the  SNdB  is  3  dB  and  DATASLCT  is  set  to  1, 
then  OATS  model  will  determine  the  angle  at  which  the  signal 
power  is  twice  the  noise  power  (3  dB).  This  is  actually 
accomplished  using  a  halving  algorithm  (7,403).  Essen¬ 
tially,  the  algorithm  establishes  the  bounds  for  the  no¬ 
offset  zenith  angle  and  cuts  the  bounds  in  half  until  it 
converges  on  the  actual  angle.  A  flag  is  set  in  the  OATS 
model  if  the  routine  does  not  converge.  This  will  result  in 
a  warning  being  printed  to  the  screen.  The  warning  does  not 
fully  define  the  problem  but  just  tells  the  modeler  that  the 
halving  routine  has  reached  40  iterations  without  converging 
to  the  answer  within  the  established  accuracy  limits.  The 
modeler  must  decide  why  the  convergence  did  not  occur. 
Frequently,  the  convergence  may  fail  because  the  angle  is 
either  larger  than  80  degrees  or  the  SNR  will  not  be 
achieved  at  any  angle. 
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Option  two  of  the  DATASLCT  field  will  calculate  the 
no-offset  zenith  angle  at  which  the  minimum  threshold  for 
detectability  is  reached.  It  uses  the  same  halving 
algorithm  as  the  SNR  computations  described  above.  However, 
the  routine  further  goes  on  to  calculate  the  time  of  track 
for  the  lidar  tracker.  This  time  of  track  algorithm  is 
built  on  multiple  assumptions.  These  are 

1)  the  satellite  is  in  a  stable  circular  orbit 

2)  the  pointing  mechanism  is  precise  enough  to  keep  the 
beam  on  the  target 

3)  there  is  no  interrupt  of  tracking  through  the 
doppler  null 

4)  there  are  no  track  interruptions  for  any  atmospheric 
anomalies  (such  as  clouds,  dust  storms  etc.) 

5)  there  are  no  interruptions  for  system  anomalies 
(such  as  noise) 

6)  the  geometry  of  track  is  symmetrical  (a  zenith  angle 
of  30  degrees  is  equivalent  to  a  zenith  angle  of  -30 
degrees ) 

7)  the  target  aircraft  remains  straight  and  level  ana 
free  from  terrain  masking 

With  these  assumptions,  option  two  returns  the  time  of 
possible  track  (i.e.,  the  signal  level  is  above  minimum 
detectable  signal  power  requirements)  in  minutes. 

An  important  note  here  is  necessary  to  show  the  major 
limitation  in  the  first  option.  Option  one  will  return  the 
zenith  angle  that  meets  the  SNdB  threshold  but  does  not 
consider  the  signal  level  at  that  angle.  The  modeler  must 
be  careful  to  check  the  power  threshold  (minimum  detectable 
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signal)  when  using  option  one  so  as  to  avoid  finding  an 
angle  where  the  SNdB  threshold  is  met  but  the  power 
threshold  is  not. 

The  final  option  in  the  data  selection  process  is 
designed  for  data  generation.  It  gives  the  modeler  the 
option  of  designating  one  of  nine  variables  as  the  variable 
to  increment  and  outputs  the  signal  power  and  signal-to- 
noise  ratio  for  each  increment  of  the  selected  variable. 

For  example,  the  modeler  could  select  to  increment  the  pulse 
power  from  1  to  10  joules  (step  size  of  1)  while  holding  all 
other  variables  constant.  For  each  incremental  value  of  the 
pulse  power,  a  data  line  is  printed.  Each  data  line 
contains  the  pulse  power,  the  receiver  power  and  the  signal- 
to-noise  ratio  for  that  incremental  run.  The  nine  options 
to  increment  are: 

1  -  no-offset  zenith  angle 

2  -  satellite  altitude 

3  -  laser  pulse  power 

4  -  laser  pulse  width 

5  -  receiver  threshold  sensitivity 

6  -  required  si gnal -to-noi se  ratio  for  detection 

7  -  diameter  of  the  receiver  collection  optics 

8  -  reflectivity  of  the  target’s  skin 

9  -  target  offset  distance  from  satellite  ground  track 

For  the  modeler’s  convenience,  the  data  output  for  this 
option  is  sent  to  a  file  called  "D.DAT”  which  is  written  to 
an  output  disk  for  later  use.  This  output  is  then  readily 
accessible  for  data  plotting  or  storage. 
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An  additional  limitation  of  the  OATS  model  is  that  it 
does  not  consider  probability  of  detection  nor  probability 
of  false  alarm.  The  data  output  and  comparisons  throughout 
the  model  are  done  on  power  thresholds  only.  The  modeler 
has  no  automated  probability  of  detection  capability. 

The  OATS  model  has  many  limitations  formed  by  the 
assumptions  made  during  its  construction.  These  limita¬ 
tions,  however,  are  soft  and  can  easily  be  altered  by  a  few 
lines  of  code.  The  core  of  the  model  is  the  lidar  range 
equation  (Equation  1)  and  is  supported  by  subfunctions  that 
can  be  manipulated  to  modify  the  program  for  additional 
analysis . 
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IV.  RESULTS 

The  basic  question  that  this  thesis  set  about  to 
answer  is  whether  or  not  a  lidar  system  on  a  satellite  could 
be  used  to  track  aircraft.  The  OATS  model  was  built  to  do 
this  analysis  with  a  wide  range  of  options  available  for 
analytical  comparisons.  If  the  OATS  model  methodology  is 
accepted  as  valid,  realistic  lidar  parameters  must  be  es¬ 
tablished  for  input  data.  For  example,  if  the  OATS  model 
were  run  using  a  pulse  power  of  10  joules  for  the  laser, 
before  the  results  can  be  claimed  as  acceptable  the  entering 
arguments  (e.g.,  10  joule  pulse)  must  be  accepted  as  valid, 
reasonable  and  realistic. 

Because  the  OATS  model  contains  so  many  variables  for 
analyzing  aircraft  trackability  from  space,  a  set  of  ac¬ 
cepted  constants  was  established.  For  the  analysis  that 
follows,  the  following  variables  are  set  up  as  constants. 


Pulsewidth 

PW 

0.1  microseconds 

Rx  Sensitivity 

RxSens 

0.001  microwatts 

Required  S/N 

SNR 

14  dB 

Processing  Gain 

PG 

3  dB 

Satel 1 i te  T rack 

SATRAK 

0  degrees 

Target  Heading 

TGTHDG 

40  degrees 

Beam  Divergence 

BMDVG 

model  #  3 

Target 

LCSMODEL 

mode  1  f  2  (brick) 

Target  size 

5  sq  meters  front 

21  sq  meters  top 
12  sq  meters  side 
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Once  the  above  variables  were  established  as  constants, 
another  set  of  variables  was  baselined  for  comparison. 

These  baseline  variables  are  listed  below  with  their  base¬ 
line  values. 


Atmosphere 

ATMOS 

model  #2  (5  km  haze) 

Zenith  Angle 

(no-offset) 

ONOZANG 

30  degrees 

Satellite  Altitude 

SATALT 

700  km 

Pulse  Power 

Joule 

1  joule 

Diameter  of  Optics 

RXD 

0.8  meters 

Tgt  Reflectance 

REFLECT 

21  * 

Tgt  Offset 

OFFSET 

100  km 

With  the  above  set  of  constants  and  baseline  values, 
the  OATS  model  was  run  for  all  four  types  of  lasers  listed 
in  the  program.  For  each  laser  type,  a  curve  plotting  the 
receiver  power  vs.  no  offset  zenith  angle  was  developed. 
These  curves  are  shown  in  Figure  12  for  the  atmosphere  with 
haze  and  in  Figure  13  for  the  atmosphere  without  haze. 

From  Figure  12,  one  can  conclude  that  holmium  yields 
the  best  return  for  the  same  system  conditions.  In  reality 
however,  many  considerations  go  into  a  determination  of 
"best"  such  as  cost,  weight,  size,  efficiencies  etc.  For 
the  sake  of  limiting  the  scope  of  this  thesis,  holmium  was 
chosen  as  the  optimum  laser  medium  for  aircraft  tracking 
strictly  on  its  performance  diagrammed  in  Figure  12.  All 
further  analysis  is  based  on  holmium,  i n  a  5  km  haze 
atmosphere . 
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RX  POWER  V8  ZENITH  ANGLE 

ISkmHAZE) 


Ce  ^  C13  o  HOIMUM  A  NEODniUM 

FIGURE  12  LASER  COMPARISON  (with  Haze) 


RX  POWER  V8  ZENITH  ANGLE 

(clear  almoepherel 


ZENITH  ANGLE 

>  C13  o  HOIMUM  A  NEODYMU 

FIGURE  13  LASER  COMPARISON  (no  Haze) 


Various  baseline  parameters  are  changed  in  the  next 
series  of  tests  to  determine  the  range  of  feasibility  for  a 
holmium  laser  in  the  5  km  atmosphere. 

The  first  parameter  to  be  varied  was  the  no  offset 
zenith  angle.  Although  this  was  accomplished  to  form  the 
holmium  curve  in  Figure  12,  Figure  14  shows  this  curve 
standing  alone.  Notice  in  Figure  14  the  "dogleg"  at  the  64 
degree  point  on  the  curve.  This  change  in  the  curve  is  a 
result  of  the  spotsize  algorithm.  From  0  degrees  to 
approximately  64  degrees,  the  spot  size  of  the  laser  was 
smaller  than  the  projected  area  of  the  target.  This  is 
reasonable  since  an  aircraft  (or  brick)  has  a  larger  proj¬ 
ected  area  from  a  top  perspective  than  a  frontal /rearward 
perspective.  After  64  degrees,  the  target’s  projected  area 
is  dominated  by  the  smaller  frontal /rearward  surface.  The 
spot  size  to  projected  area  ratio,  x.  becomes  one  and  no 
longer  limits  the  LCS  equations  in  the  model.  The  CATS 
model  lists  the  value  for  x  as  the  program  is  executing  and 
therefore  gives  the  modeler  an  insight  on  the  occurrence  of 
"doglegs'"  prior  to  plotting  data. 

During  the  next  set  of  OATS  model  runs,  the  satellite 
altitude  was  allowed  to  vary.  By  cueing  on  available 
tracking  time  (TRAKTIME),  the  optimum  satellite  altitude  was 
computed.  This  is  graphically  shown  in  Figure  15.  Recall 
that  track  time  is  the  time  during  a  single  satellite  orbit 
that  the  target  is  trackable  (within  minimum  detectable 
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signal  thresholds).  For  the  baseline  that  was  previously 
I  established,  the  maximum  available  single  orbit  track  time 

is  approximately  23  minutes  for  a  satellite  at  approximately 
3500  km  altitude. 

m  The  plot  of  satellite  altitude  vs.  receiver  power  is 

shown  in  Figure  16.  As  is  predictable,  the  receiver  power 
drops  off  as  the  distance  increases.  Again,  however,  one 
can  see  an  apparent  anomaly  at  approximately  1600  km  on  the 
plot.  This  discontinuity  in  the  curve  is  attributable  again 
to  the  spotsize  correction  algorithm.  As  the  spot  size 
*  grows  beyond  the  projected  area  of  the  target,  the  power 

received  at  the  satellite  changes  at  a  different  rate  also. 
Beyond  approximately  1600  km,  the  beam  becomes  a  'flood 
B  light",  totally  illuminating  the  target’s  silhouette. 

Before  1600  km,  the  beam  was  smaller  that  the  target’s 
silhouette  and,  therefore,  only  illuminated  a  patch"  on  the 
®  target. 

The  next  set  of  data  shows  the  change  in  receiver  power 
as  the  offset  distance  increases.  Recall  that  the  offset 
distance  is  the  distance  from  the  satellite’s  ground  track 
to  the  target,  always  a  perpendicular  projection  from  the 
satellite’s  ground  track  line.  This  is  the  smallest  dis¬ 
tance  from  the  target  to  the  satellite’s  ground  track.  From 
the  plot,  Figure  17,  one  can  see  that  the  change  appears 
linear.  The  relationship  is  not  linear  since  the  vertical 
'  axis  is  logarithmic.  Establishing  a  minimum  detection  level 
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OFFSET  (km) 

FIGURE  17  Rx  POWER  V8  OFFSET 


at  0.001  microwatts  (-60  dBm  on  the  scale),  it  appears  that 
the  lidar  tracking  system  could  easily  track  targets  of  this 
size  out  as  far  as  1000  km  from  the  satellite’s  ground  po¬ 
sition.  This  does,  however,  assume  a  very  clean  atmosphere. 

The  next  parameter  of  variation  is  target  reflectance. 
Measured  as  a  percentage  of  power  reflected,  this  parameter 
is  a  function  of  the  target’s  skin,  whether  paint  or  bare 
skin  material.  This  is  graphically  displayed  in  Figure  18. 
Notice  that  the  rate  of  change  of  the  curve  in  the  lower 
reflectance  region  (under  20%  reflectance)  is  very  large 
while  the  rate  of  change  beyond  50%  is  relatively  constant. 
This  curve  reveals  two  significant  points  to  consider  if 
such  a  lidar  tracking  system  were  to  be  built.  The  first  is 
that  targets  that  purposely  wanted  to  be  tracked  should  use 
a  paint  that  is  highly  reflective  at  the  wavelength  in 
question.  If  a  paint  with  a  reflectance  value  of  80%  were 
available,  the  target  could  have  a  receiver  power  30  dB 
above  the  minimum  detectable  signal  level  (using  these 
optimistic  baseline  values).  This  produces  a  large  track- 
ability  domain  for  the  target.  The  other  point  to  consider 
is  that  targets  that  do  not  want  to  be  tracked,  could 
greatly  complicate  the  tracking  problem  by  using  a  highly 
absorptive  paint.  Each  small  increase  in  absorptivity  below 
10%  yields  a  great  decrease  in  the  amount  of  reflected 
power.  Note,  however,  that  even  at  a  reflectivity  of  0.04%, 
the  defined  target  is  still  trackable  at  a  no-offset  zenith 
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angle  of  30  degrees  and  an  offset  of  100  km.  This  shows 
that  the  window  of  trackabi lity  is  still  relatively  large. 

A  30  degree  no  offset  zenith  angle  results  in  a  track  time 
of  approximately  1.8  minutes,  more  than  enough  to  establisn 
a  target’s  position,  heading,  and  speed. 

The  next  two  variables  that  were  changed  for  analysis 
are  functions  of  the  lidar  system  itself.  These  parameters 
are  the  laser  pulse  power  or  pulse  energy  and  the  diameter 
of  the  optics. 

The  pulse  energy  was  varied  over  two  orders  of 
magnitude.  This  range  shows  a  very  significant  result  in 
the  under-2  joules  region.  In  this  region,  small  changes  in 
pulse  energy  resulted  in  dramatic  change  in  receiver  power. 
Above  4  joules,  the  small  increases  in  pulse  power  yield 
relatively  smaller  increases  in  receiver  power.  Again, 
notice  that  even  at  small  pulse  energy  values,  the  track- 
ability  of  the  target  is  large.  The  receiver  power  is  more 
than  10  dB  above  the  minimum  detectable  threshold  for  pulse 
energy  values  as  low  as  0.1  joules.  The  receiver  power  is 
plotted  against  the  pulse  energy  in  Figure  19. 

Figure  20  plots  the  diameter  of  the  optics,  whether 
mirrors  or  lenses,  against  the  receiver  power.  As  opposed 
to  the  other  plots,  this  plot  shows  a  very  large  change  in 
receiver  power  for  small  changes  in  diameter  of  the  optical 
system.  The  plot  shows  that  certain  diameters  yield  re¬ 
ceiver  powers  below  the  minimum  detectable  threshold  (set  to 
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-60  dBm).  The  curve  levels  off  after  about  1  meter  diameter 
with  a  slope  showing  approximately  3  dB  of  gain  for  each  50 
cm  in  the  diameter.  Small  diameter  optics  do  not  appear  to 
be  practical  since  the  receiver  power  drops  rapidly  under 
0.5  meters  diameter.  For  spaceborne  lidar  systems  however, 
size  and  weight  are  extremely  important  and  very  large 
diameter  systems  are  also  constrained  by  the  space  vehicle 
constraints . 
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CONCLUSIONS  AND  RECOMMENDATIONS 

From  the  results,  one  could  hastily  conclude  that 
holmium  lidar  tracking  of  aircraft  is  easy  and  applicable 
over  a  wide  range  of  system  and  target  characteristics. 
However,  the  results  must  be  viewed  from  a  proper  perspec¬ 
tive.  The  assumptions  upon  which  the  OATS  model  is  built 
are  crucial  and  directly  determine  the  results  of  the  model. 
The  assumptions  must  therefore  always  be  included  in  any 
discussion  of  the  results. 

In  general,  the  assumptions  which  were  made  in 
constructing  the  OATS  model  are  "optimistic",  resulting  in 
better  detectivity  than  can  be  expected  in  actual  use. 
However,  even  with  the  optimism,  this  study  does  show  that 
under  relatively  ideal  conditions,  lidar  tracking  of 
aircraft  from  space  may  be  feasible.  Further  study  is 
required  to  fully  develop  the  feasibility  envelope  for  such 
a  system. 

The  areas  of  further  study  in  this  field  can  be 
clustered  into  three  main  categories.  These  areas  are: 

1  -the  satellite  system  itself,  which  includes  such 

areas  as  the  laser,  receiver,  orbital  control, 

tracking  and  processing  algorithms,  etc. 

2  -the  environment,  which  includes  further  atmospheric 

attenuation  analysis  for  more  types  of  atmospheres 

3  -the  target,  which  includes  further  study  of  aircraft 

reflection  characteristics. 


Some  specific  area  of  improvement  are  therefore 
presented  for  future  study  and/or  improvement  for  the  OATS 
model . 

The  first  recommended  improvement  to  the  OATS  model  is 
to  obtain  FASCODE  atmospheric  attenuation  coefficients 
commensurate  with  an  atmosphere  containing  some  clouds. 

This  improvement  could  be  used  to  estimate  how  often  a 
target  could  be  expected  to  be  trackable.  For  now,  the 
results  show  excellent  results  for  cooperative  atmospheres 
only,  the  noncooperative  atmospheres  need  to  be  examined. 

The  next  improvement,  alluded  to  ea’'lier,  is  to  correct 
the  ranging  algorithm  for  optical  path  length  instead  of  the 
direct  path  length  assumption  used  within.  The  maximum 
error  for  this  approximation  in  the  far  IR  wavelengths  is 
approximately  2%  (0  =  90*’).  The  quickest  method  to  achieve 
better  ranging  would  be  to  use  the  range  value  already 
established  in  FASCODE.  These  values,  taken  when  the 
attenuation  values  are  taken,  can  be  fit  to  a  curve  (just  as 
the  transmittance  values  are)  and  incorporated  into  the  OATS 
model  as  a  fifth-order  polynomial  (just  like  transmittance). 

The  next  improvement  is  the  addition  of  doppler  con¬ 
siderations.  The  structure  in  the  OATS  model  already  estab¬ 
lishes  the  geometry  for  doppler  calculations.  A  window 
defining  doppler  tracking  requirements  (i.e.  amount  of 
doppler  shift  necessary  for  tracking)  could  be  added  to  show 
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where  the  doppler  null  would  preclude  tracking  and  define 
where  tracking  would  be  re-established. 

To  maintain  the  realism  of  the  model,  further  noise 
sources  should  be  added.  In  doing  this,  one  may  decide  to 
include  actual  receiver  configurations,  materials  and 
circuitry  so  as  to  get  a  better  noise  figure.  Different 
configurations  could  be  added  as  options  into  the  model  just 
as  different  lasers  at  specific  resonances  are  currently 

/ 

i ncorporated . 

Another  consideration  that  should  be  considered  for 
further  work  in  this  area  is  to  re-establish  the  coordinate 
system  for  the  geometry  to  a  geocentric  reference  frame. 

This  reference  frame  would  provide  an  opportunity  to  model 
several  satellites,  several  targets  and  do  more  complex 
multi-target  and  multi-satellite  analysis.  More  exact 
analysis  could  be  performed  with  offset  targets.  Currently, 
the  OATS  model  is  limited  in  its  accuracy  as  the  offset 
distances  grow  larger.  Precise  spherical  geometry  could  be 
used  to  compute  accurate  geometries  for  large  offset 
distances . 

One  of  the  possible,  albeit  difficult  to  obtain, 
improvements  to  the  OATS  model  would  be  to  add  actual 
aircraft  lidar  cross  section  data  along  with  actual  skin 
paint  reflectance  values.  If  this  data  were  available,  the 
applicability  of  the  model  to  actual  aircraft  would  be 
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greatly  enhanced.  It  would  also  tend  to  validate  or  refute 
the  flying  brick  model.  It  may  even  reveal  a  potential 
improvement  for  the  flying  brick  model  that  would  increase 
its  applicability. 

The  next  bundle  of  recommendations  relate  to  side 
investigations  and  are  not  necessarily  improvements  to  OATS. 
The  first  of  these  investigations  would  be  the  study  of  the 
potential  use  of  this  lidar  system  in  a  bistatic  scenario. 
Bistatic  geometries  could  add  another  dimension  to  the  sys¬ 
tem  capabilities  that  were  lost  from  clouds  or  non-diffuse 
reflections. 

This  proposed  orbital  lidar  system  would  need  very  fine 
beam  steering  capabilities  to  accomplish  its  assigned  tasks. 
Analysis  of  precision  high  power  beam  steering  could  be  the 
"green  light"  or  the  "red  light"  for  an  entire  program  like 
an  orbital  lidar.  If  no  technology  exists  to  steer  the  beam 
in  the  directions  required,  then  the  system  is  impractical. 
Along  with  beam  steering,  one  would  need  to  develop  scan 
algorithms  for  searching,  acquiring  and/or  tracking.  A 
simulation  model  would  be  an  excellent  search,  acquisition 
and  track  evaluation  tool  and  could  use  the  OATS  model  as  a 
subroutine  to  the  main  simulation  program. 

Another  potential  simulation  for  system  design  would  be 
to  design  a  model  to  test  different  satellite  constellations 
for  multiple  orbital  lidars.  One  would  establish  a  certain 
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criterion  for  coverage  and  proceed  to  build  up  constellation 
configurations  to  meet  the  coverage  requirements.  Again, 
this  one  study  could  show  the  tracking  from  space  concept  as 
infeasible  due  to  sheer  number  of  satellites  required.  Then 
again,  it  may  show  that  orbital  lidar  systems  are  maximized 
for  certain  satellite  conf igurations. 

The  final  and  least  precise  study  that  must  be  accom¬ 
plished  is  a  cost  analysis.  This  would  follow  the  config¬ 
uration  study.  Once  the  configuration  study  revealed  the 
number  of  satellites  required,  then  a  cost  for  a  space-based 
tracking  system  could  be  assessed. 

A  final  recommendation  here  is  from  an  idea  that  Dr. 
Paul  McManamon  had  during  a  discussion.  A  space-based  lidar 
system  that  was  built  as  a  sensor  for  the  Strategic  Defense 
Initiative  could  potentially  be  used  for  tracking  aircraft 
as  well  (18).  Expanding  on  the  multi-use  idea,  one  can 
easily  find  scores  of  uses  for  a  lidar  in  space.  Once  in  an 
established  orbit,  it  could  be  used  to  track  shuttle 
missions,  track  satellites,  and  even  space  debris.  The  same 
system  could  be  used  for  atmospheric  and  environmental 
study.  For  example,  an-orbital  lidar  could  potentially 
track  migrating  Canadian  geese  and  use  the  data  for  aircraft 
avoidance  or  biological  studies. 


APPENDIX  A  OATS  MODa  CODE 


’  OPTICAL  AIRCRAFT  TRACKING  SYSTEM  (OATS)  MODEL  *»*» 

•  —  _ I _ , _ 

DECLARE  SUB  MENU  ( ) 

DECLARE  SUB  CHGMENU  ( ) 

DECLARE  SUB  CALC  () 

DECLARE  FUNCTION  BETA!  () 

DECLARE  FUNCTION  dBm!  (PWX) 

DECLARE  FUNCTION  LCS  (ZANG) 

DECLARE  FUNCTION  TRANSMIT!  () 

DECLARE  FUNCTION  RANGE!  (SAT ALT) 

DECLARE  FUNCTION  SPOTS2!  (BMOVG.  RNG) 

DECLARE  FUNCTION  TRAKTIME!  () 

DECLARE  FUNCTION  TPaREA!  () 

COMMON  SHARED  ATTEN.  ATMOS,  N02ANG,  20RNG,  R,  Al,  A3,  A2,  XI,  OFFSET,  SATALT. 
— ZANG 

COMMON  SHARED  JOULES,  ONOZANG,  T,  RNG,  TLCS,  SZ,  BMOVG,  RXAREA,  LCSMODEL, 

— RADSPHR 

COMMON  SHARED  CIU,  C2tl,  C3«,  C4»,  C5»,  CI«,  PW.  WL,  PWX,  S,  PG,  Z,  LASER. 

— BEAMOO 

COMMON  SHARED  SNdB,  OATASLCT,  SATRAK,  TGTHOG,  PWI,  RXmlcro,  PGdB,  RXD,  REFLECT 
COMMON  SHARED  SLCTR,  RA1,  RA2,  RA3,  RXSENS,  SNR,  CHKR1,  CHKR2,  CHI 
OPEN  "O.OAT"  FOR  OUTPUT  AS  #1 

4.  , 

COLOR  2 

PRINT  "  WELCOME  TO  THE  OPTICAL  AIRCRAFT  TRACKING  SYSTEM  MODEL  (OATS)  " 

PRINT 

PRINT  " - VERSION  1.1  21  NOV  1988 - ' 

PRINT 

PRINT  "OATS  IS  a  comouter  model  of  a  laser  radar  (LIDAR)  system  used  to  track 
PRINT  "aircraft.  It  offers  a  wide  range  of  ootions  for  engineering  and 
PRINT  "scientific  use  with  a  variety  of  output  options.  Specific  assumptions' 
PRINT  "imbedded  in  this  model  are  explained  in  the  thesis  document  that  goes" 
PRINT  "with  this  program. 

COLOR  3 


PRINT 

PRINT 

PRINT  " - CAUTION - 

PRINT  "  Use  of  this  program  without  understanding  the  assumptions" 
PRINT  ”  upon  which  it  is  built  may  cause  undo  embarassment  &  harm" 
PRINT  "  YOU  WERE  WARNED  !" 

PRINT  - 


COLOR  2 
PRINT 
PRINT 
PRINT 


PRINT  "  This  version  of  OATS  Is  specifically  designed  for  space-based 
—applications" 

PRINT  "  It  assumes  a  target  at  low  altitude  (10  meters)  and  offers  four 
—different 

PRINT  "  laser  types  to  chose  from.“ 

LOCATE  25,  21 

INPUT  "PRESS  RETURN  TO  CONTINUE";  PR$ 

CLS 

INPUT  “ENTER  SATELLITE  ALTITUDE  IN  KILOMETERS";  SATALT 
SELECT  CASE  SATALT 


CASE  IS  <  120 
CLS 

PRINT  " - WARNING  WARNING  WARNING - 

PRINT  '• -  THIS  IS  PRETTY  LOW  -  " 

PRINT 

CASE  IS  >  30000 
CLS 

PRINT  ■ - WARNING  WARNING 

— WARNING - " 

PRINT  " -  THIS  IS  VERY  HIGH  FOR  A  TRACKING  SATELLITE 


PRINT 
CASE  ELSE 
PRINT 

END  SELECT 
CHOSELASER:  PRINT 

PRINT  "ENTER  THE  NUMBER  FOR  THE  TYPE  OF  LASER  THAT  YOU  DESIRE  TO  TEST" 
PRINT  "  1  -  C02  with  normal  Carbon  12  molecules' 

PRINT  "  2  -  C02  with  Carbon  13  (Isotope)  molecules" 

PRINT  "  3  -  HOLMIUM/YLF  centered  at  2.06  micrometers" 

PRINT  "  4  -  NEOOYNIUM  in  silicate  glass'” 


"Wavelength  data  from  the  following  sources 

’C02  wavelengths  from  "The  C02  Laser"  by  W.J.Witteman 

’  Berl inzSpringer-Verlag,  1987 

’  P20  line  for  C-12 - page  24 

’  P20  line  for  C-i3 - page  27 

'Holmium  wavelength  from  Ken  Schepler,  AFWAL 
’  telephone  conversation  13  Oct  1988 

'Neodymium  wavelength  from  "CRC  Handbook  of  Laser  Science  & 

’  Technology  Volume  V:  Optical 

’  Materials  Part  3:  Applications, 

’  Coatings  and  Fabrication" 

’  Editor  Marvin  J.  Weber 

’  Boca  Raton  Florida: CRC  Press  Inc. 

’  1987  — page  363, 

’  TYPE-  Silicate  glass 

’  Glass  Designation  -  LG-660 

’  Manufacturer  -  Schott 

INPUT  LASER 
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SELECT  CASE  LASER 
CASE  1 

WL  =  10.59104  /  1000000 
PRINT 

CASE  2 

WL  =  11.1494  /  1000000 
PRINT 

CASE  3 

WL  =  2.06  /  1000000 
PRINT 

CASE  4 

WL  =  1.057  /  1000000 
PRINT 
CASE  ELSE 
PRINT 

PRINT  "INPUT  ERROR  - TRY  AGAIN  WITH  A  GOOD  VALUE" 

GOTO  CHOSELASER 
ENO  SELECT 
PRINT 

INPUT  "ENTER  LASER  PULSE  ENERGY  (In  Joules)";  JOULES 
PRINT 

INPUT  "ENTER  PULSEWIOTH  OF  EACH  LASER  PULSE  (in  microseconds)";  PWI 

PW  =  PWI  /  1000000 

PRINT 

PRINT 

PRINT  "NOTE:  Model  assumes  a  Bandwidth-pul  sew idth  product  of  1  for  S/N 
—calculations" 

PRINT 


PRINT 

INPUT  "ENTER  REQUIRED  S/N  VALUEdn  dB)";  SNdB 
SNR  =  10  '  (SNdB  /  10) 

PRINT 

INPUT  "ENTER  THE  RECEIVER  SENSITIVITY  (in  Microwatts)";  RXmicro 

RXSENS  =  RXmicro  /  1000000 

PRINT 

INPUT  "ENTER  PROCESSING  GAIN  (in  dB)";  PGdB 
PG  =  10  "  (PGdB  /  10) 

PRINT 

INPUT  "ENTER  EFFECTIVE  DIAMETER  OF  THE  RECEIVING  OPTICS  (in  meters)  ';  RXO 

RXAREA  =  (3.141592651  *  (RXO  *  2))  /  4 

OIVG: 

PRINT 

PRINT  "ENTER  LASER  BEAM  MODEL  FOR  THE  BEAM  DIVERGENCE" 

PRINT  "  1  -  SPECIFY  3dB  BEAM  DIVERGENCE  MANUALLY" 

PRINT  "  2  -  USE  RCVR  DIAMETER  &  3dB  DIFFRACTION  LIMIT" 

PRINT  "  3  -  USE  RCVR  DIAMETER  &  WORST  THAN  DIFFRACTION  LIMIT" 

INPUT  BEAMOD 

SELECT  CASE  BEAMOD 
CASE  1 


INPUT  ENTER  BEAM  DIVERGENCE  (In  Radians)";  BHOVG 
CASE  2 

BHOVG  =  1.08  «  WL  /  RXO 
CASE  3 

BHOVG  =  1.3  «  WL  /  RXO 
CASE  ELSE 

PRINT  "TRY  AGAIN,  INPUT  VALUE  WAS  IN  ERROR  ' 

GOTO  DIVG 
END  SELECT 

GaInTX  =  12.56637062*  /  (BHOVG  "  2) 

CHOSEATHOS:  PRINT  " 

PRINT 

PRINT  “ENTER  NUHBER  FOR  THE  ATHOSPHERIC  HOOEL  DESIRED" 

PRINT  “  1  -  CLEAR,  no  aerosols" 

PRINT  "  2  -  HAZE,  5  km  visibility  rural  type  aerosol" 

PRINT  "  3  -  SELECTED  ATTENUATION,  manually  entered  attenuation" 

INPUT  ATHOS 

SELECT  CASE  ATHOS 
CASE  1,  2 
PRINT 
CASE  3 

PRINT  "  ENTER  THE  ATTENUATION  TO  ADD  TO  THE  ALREADY  EXISTING" 

PRINT  "  ATTENUATION  FROH  THE  HOLECULAR  ABSORPTION  (in  dB)" 

PRINT  "  i.e.  50  X  attenuation  from  an  aerosol  would  be  3dB" 

PRINT  "  entry  is  for  ONE-WAY  passage  through  the  atmosphere" 

PRINT  "  at  a  zenith  angle  of  ZERO  (straight  up  through  the 
—atmosphere)" 

PRINT  "  — an  exponential  curve  (trans  vs  range)  is  used  for  all" 
PRINT  "  zenith  angles  other  than  zero." 

INPUT  ATTENdB 

ATTEN  =  10  *  (ATTENdB  /  10) 

CASE  ELSE 
PRINT 

PRINT  "INPUT  ERROR  - TRY  AGAIN  WITH  A  GOOD  VALUE 

GOTO  CHOSEATHOS 

END  SELECT 
CHOSEHODEL:  PRINT 

PRINT  "ENTER  TARGET  CROSS  SECTION  HODEL" 

PRINT  "  1  -  TARGET  IS  A  SPHERE" 

PRINT  "  2  -  TARGET  IS  A  BRICK" 

INPUT  LCSHODEL 

SELECT  CASE  LCSHODEL 
CASE  1 

INPUT  "WHAT  IS  THE  RADIUS  OF  THE  SPHERE  IN  METERS":  RADSPHR 
PRINT  "WHAT  IS  THE  REFLECTANCE  OF  THE  SPHERE  ?" 

PRINT  "reflectance  is  that  percentage  of  energy  at  the 
PRINT  "wavelength  of  interest  that  is  reflected  (between  0  TO 

— 100)" 

INPUT  REFLECT 


R  =  REFLECT  /  100 
CASE  2 

PRINT  "THE  FLYING  BRICK  MODEL  IS  A  THREE  DIMENSIONAL  MODEL" 

PRINT  "OF  AN  OBJECT  SIMPLIFYING  THE  TARGET  TO  THREE  ORTHOGONAL 
PRINT  "PLANES  LIKE  THREE  SIDES  OF  A  BRICK.  THEREFORE,  THE  THREE" 
PRINT  "FACES  OF  THE  SYMETRICAL  BRICK  REPRESENT  THE  PROJECTED 

— AREA" 

PRINT  "OF  THE  TARGET  FROM  A  FRONTAL/REAR  VIEW,  FROM  A  TOP  VIEW 

— AND" 

PRINT  "FROM  A  SIDE  VIEW  ...  essentially  the  target's  silhouette  " 
INPUT  "ENTER  THE  FRONT  /  REAR  PROJECTED  AREA  (sqr  meters)";  RA1 
INPUT  "ENTER  THE  TOP  PROJECTED  AREA  (sqr  meters)';  RA2 
INPUT  "ENTER  THE  SIDE  PROJECTED  AREA  (sqr  meters)";  RA3 
PRINT  "WHAT  IS  THE  REFLECTANCE  OF  THE  SKIN  OF  THE  TARGET  ?" 

PRINT  "reflectance  is  that  percentage  of  energy  at  the  " 

PRINT  "wavelength  of  interest  that  is  reflected  (between  0  to 

— 100)" 

INPUT  REFLECT 

R  =  REFLECT  /  100 
CASE  ELSE 

PRINT  »«««««««««««««««««««»»*««««" 

PRINT 

PRINT  "  -  INPUT  MUST  BE  1  or  2 . TRY  AGAIN  - " 

PRINT 

PRINT 

GOTO  CHOSEMODEL 
END  SELECT 


HDGINPUT:  PRINT 

INPUT  "ENTER  SATELLITE  GROUND  TRACK  HEADING  REFERENCING  TRUE  NORTH  (degrees)"; 
SATRAK 

INPUT  "ENTER  TARGET  HEADING  REFERNCING  TRUE  NORTH  (degrees)";  TGTHDG 
XXXX  =  BETA 

IF  CHKR1  =  1  GOTO  HDGINPUT 
IF  CHKR2  =  1  GOTO  HDGINPUT 

PRINT  "The  offset  distance  is  the  smallest  distance  from  the  target  to  the 
— satellite’s" 

PRINT  "ground  track.  It  therefore  scribes  a  line  that  is  orthogonal  with  the' 
PRINT  "satellites  ground  track" 

INPUT  "ENTER  OFFSET  DISTANCE  (in  kilometers)";  OFFSET 
PRINT 

PRINT  "ENTER  TYPE  OF  DATA  DESIRED" 

PRINT  "  1  -  Calculate  the  NOZANG  where  the  S/N  meets  required  level" 

PRINT  "  2  -  Calculate  the  NOZANG  where  the  minimum  detectible  pwr  is  met" 

PRINT  "  3  -  Make  one  parameter  vary  and  calculate  min  pwr  &  S/N  over  a 

—range" 

INPUT  DATASLCT 
ALPHA: 


A1  =  ABS(RA1  *  (COS(BETA))  '  2)  +  (RA3  *  (SIN(BETA))  ‘  2) 
A2  =  RA2 

A3  =  ABS(RA3  *  (COS(BETA))  “  2)  +  (RAl  *  (SIN(BETA))  *  2) 


SELECT  CASE  OATASLCT 
CASE  1 

RIGHT  =  80:  LEFT  =  0:  MIO  =  40 
COUNT  =  0 

LAP1:  COUNT  =  COUNT  ♦  1:  ONOZANG  =  MIO:  NOZANG  =  DNOZANG  /  57.2957795m 
IF  COUNT  =  40  THEN 
CLS 

PRINT  " - 

PRINT  "  ***  WARNING  »**  " 

PRINT  "  count  has  reached  40  without  converging  to  an  answer 
PRINT 


GOTO  QUIT1 
END  IF 
CALC 

delta  :  S  -  SNR 
SELECT  CASE  delta 
CASE  IS  >  .0001 
LEFT  s  MID 

MID  =  (LEFT  +  RIGHT)  /  2 
GOTO  LAP1 
CASE  IS  <  -.0001 
RIGHT  =  MIO 

MIO  =  (LEFT  +  RIGHT)  /  2 
GOTO  LAP1 
CASE  ELSE 

PRINT  "+++++++♦++++++++++  BINGO  YOGI  ++++++++++++++-f++++++++++ 
END  SELECT 
QUIT1:  PRINT 
PRINT 

PRINT 

PRINT  "NO  OFFSET  ZENITH  ANGLE  FOR  THE  REQUIRED  SIGNAL  TO  NOISE  RATIO  IS" 

PRINT  ■■  ",  MIO 

PRINT 


RIGHT  =  80:  LEFT  =  0:  MID  =  40 
COUNT  =  0 


LAP2:  COUNT  =  COUNT  +  1:  ONOZANG  =  MID:  NOZANG  =  ONOZANG  /  57.29577951* 
IF  COUNT  r  40  THEN 
CLS 


- - - -  - 

PRINT  ■■  ***  WARNING  ***  ' 

PRINT  "  count  has  reached  40  without  converging  to  an  answer 
PRINT 

GOTO  QUIT2 
£N0  IF 
CALC 

delta  =  PWX  -  RXSENS 

PRINT  «1,  USING  "§».»»»» - ;  DNOZANG,  N02ANG,  PWX,  RXSENS,  delta 

PRINT  delta 

SELECT  CASE  delta 

CASE  IS  >  .00000000000001* 

LEFT  =  MID 

MID  =  (LEFT  +  RIGHT)  /  2 
GOTO  LAP2 

CASE  IS  <  -.00000000000001* 

RIGHT  =  MID 

MID  =  (LEFT  +  RIGHT)  /  2 
GOTO  LAP2 
CASE  ELSE 
CLS 
PRINT 

PRINT  ■++++++++++++++++++  8INGO  -  CONVERGENCE  COMPLETE  ! 
~+++++++++++++++++++++++++  ■■ 

END  SELECT 
QUIT2:  PRINT 

IF  S  >  SNR  THEN 
PRINT 

PRINT 

PRINT  "NO  OFFSET  ZENITH  ANGLE  WHERE  SIGNAL  POWER  IS  OETECTIBLE  IS" 
PRINT  “  ",  MID,  "Degrees" 

PRINT 

PRINT 

PRINT  "  The  track  time  available  with  this  zenith  angle  and  satellite 
— altitude" 

PRINT  "  ",  TRAKTIME,  "Minutes" 

PRINT 

ELsi 

CLS 

PRINT 

PRINT 

PRINT  "SIGNAL  TO  NOISE  RATIO  IS  INSUFFICIENT  FOR  TARGET  TRACKING" 

END  IF 


PRINT  THE  VARIABLES  THAT  YOU  ARE  ALLOWED  TO  VARY  ARE" 

PRINT  "  1  -  ONOZANG  (th«  No  Offset  Zenith  Angle  In  Degrees)" 

PRINT  "  2  -  SATALT  (SATellite  ALTitude  In  k1 lometers)" 

PRINT  "  3  -  JOULES  (the  laser  pulse  power  in  JOULES)" 

PRINT  "  4  -  PW  (the  Pulse  Width  In  microseconds)" 

PRINT  "  5  -  RXSENS  (RX  minimun  detectible  signal  level  In 

--microwatts)" 

PRINT  "  6  -  SNR  (the  required  S/N  to  declare  a  target)" 

PRINT  "  7  -  RXD  (the  Diameter  of  the  Receiver  collection  optics)" 

PRINT  “  8  -  REFLECT  (the  REFLECTIvlty  of  the  target’s  skin  In 

—percent)" 

PRINT  "  9  -  OFFSET  (distance  from  satellite’s  ground  track  to 

— target ) ’ " 

PRINT 

INPUT  "ENTER  THE  OPTION  YOU  WISH  TO  VARY";  DCHOICE 
PRINT 

INPUT  "ENTER  LOW  BOUND  FOP  YOUR  CHOICE";  BOTTOM 
PRINT 

INPUT  "ENTER  UPPER  BOUND  FOR  YOUR  CHOICE";  TOP 
PRINT 

INPUT  "ENTER  STEP  SIZE";  INCR 
SELECT  CASE  DCHOICE 

CASE  2,  3,  4,  5,  6,  7,  8,  9 

INPUT  "AT  WHAT  NO  OFFSET  ZENITH  ANGLE  00  YOU  DESIRE  THE  DATA 
—(degrees)";  ONOZANG 
CASE  ELSE 
PRINT 
END  SELECT 
SELECT  CASE  DCHOICE 
CASE  1 

FOR  ONOZANG  =  BOTTOM  TO  TOP  STEP  INCR 
CALC 

PRINT  *1,  USING  - ;  ONOZANG,  PWX,  dBm(PWA).  S 

PRINT  USING  - ;  ONOZANG,  RNG,  PWX,  dBin(PWX),  S 

NEXT 
CASE  2 

FOR  SATALT  =  BOTTOM  TO  TOP  STEP  INCR 
CALC 

PRINT  USING  "$$.$$$» - ;  SATALT,  PWX,  dBm(PWX),  S 

PRINT  «1,  USING  "»$.»$$$ - ;  SATALT,  PWX,  dBm(PWX),  S 

NEXT 
CASE  3 

FOR  JOULES  =  BOTTOM  TO  TOP  STEP  INCR 
CALC 

PRINT  USING  "«*.•»»* - ;  JOULES,  PWX,  dBm(PWX),  S 

PRINT  #1,  USING  "**.***» - ;  JOULES,  PWX,  dBm(PWX),  S 

NEXT 
CASE  4 

FOR  PWI  :  BOTTOM  TO  TOP  STEP  INCR 

PW  :  PWI  /  1000000 

CALC 
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PRINT  USING  - ;  PWl,  PWX,  dBm(PWX),  S 

PRINT  #1,  USING  . ;  PWI,  PWX,  dBmCPWX),  S 

NEXT 
CASE  5 

FOR  RXmlcro  =  BOTTOM  TO  TOP  STEP  INCR 

RXSENS  =  RXfflicro  /  1000000 

CALC 

PRINT  USING  - ;  RXmicro,  PMX,  dBm(PWX),  S 

PRINT  i1,  USING  “•«.#•«# - ;  RXimcro,  PWX,  dBin(PWX),  S 

NEXT 
CASE  6 

FOR  SNdB  =  BOTTOM  TO  TOP  STEP  INCR 
SNR  =  10  ‘  (SNdB  /  10) 

CALC 

PRINT  USING  . ;  SNdB,  PWX,  dBm(PWX),  S 

PRINT  i1,  USING  - ;  SNdB,  PWX,  dBm(PWX).  S 

NEXT 
CASE  7 

FOR  RXO  =  BOTTOM  TO  TOP  STEP  INCR 
RXAREA  s  (3.14159265i  »  (RXD  "  2))  /  4 
SELECT  CASE  BEAMOO 
CASE  2 

BMDVG  :  1.08  *  WL  /  RXO 
CASE  3 

BMOVG  :  1.3  «  WL  /  RXO 
CASE  ELSE 

PRINT  "CHECK  THE  VARIATION  IN  YOUR  BEAM  DIVERGENCE  AS  RXD  CHANGES" 
ENO  SELECT 
CALC 

PRINT  USING  "$*.»»»$ - ;  RXD,  PWX,  dBm(PWX),  S 

PRINT  *1,  USING  - ;  RXO,  CHI,  dBn>(PWX),  S 

NEXT 
CASE  8 

FOR  REFLECT  =  BOTTOM  TO  TOP  STEP  INCR 

R  =  REFLECT  /  100 

CALC 

PRINT  USING  "Hi.*##*"'  REFLECT,  PWX,  dBm(PWX),  S 

PRINT  *1,  USING  ■’*«.*»«# - ;  REFLECT,  PWX,  dBm(PWX),  S 

NEXT 
CASE  9 

FOR  OFFSET  =  BOTTOM  TO  TOP  STEP  INCR 
CALC 

PRINT  USING  •■*».#««* . ;  OFFSET.  PWX,  dBin(PWX),  S 

PRINT  *1,  USING  - ;  OFFSET,  PWX,  dBm(PWX),  S 

NEXT 

CASE  ELSE 

PRINT  "CHOICE  WAS  IMPROPERLY  ENTERED. .. DATA  IS  QUESTIONABLE" 

ENO  SELECT 
CASE  ELSE 

PRINT  "SELECTION  OF  DATA  WAS  IMPROPERLY  ENTERED. .. DATA  IS  QUESTIONABLE" 
ENO  SELECT 
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COLOR  17,  7 
LOCATE  25,  12 

INPUT  "  +■♦>++++++++++++  PRESS  RETURN  TO  CONTINUE  ++++++■*■•*■ 

COLOR  2,  0 

HENU 

IF  SLCTR  =  0  GOTO  ALPHA 
CLS 

LOCATE  14,  12 
COLOR  0,  2 

PRINT  "THAT  COMPLETES  THE  OATS  MODEL  FOR  NON - i.e.  THE  END" 


+•*■++++-*■•*■+■*•+++'■;  PR$ 


FUNCTION  BETA 

SELECT  CASE  SATRAK 
CASE  IS  <  0 
CLS 

LOCATE  12,  1 

PRINT  ''«<««*«*««*«***«**««*»« «*t <«*<««« 

PRINT  "ERROR  ERROR  ERROR  -  DATA  IS  OUT 

PRINT  " - Re-enter  satellite  track-- 

PRINT  " - Re-enter  target  heading— 

PRINT  '’*****t*t**t»**t*ti**tttt*tttt**tM**t*t*t** 
CHKR1  =  1 
CASE  IS  <:  180 
CHKR1  :  0 
CASE  IS  >  360 
CLS 

LOCATE  12,  1 

PRINT  •■**»*»* t«*«**»»»»*s»**» »«»»»**»»»*«»*«*»»$ 

PRINT  "ERROR  ERROR  ERROR  -  DATA  IS  OUT 

PRINT  " - Re-enter  satellite  track — 

PRINT  " - Re-enter  target  heading-- 

PRINT  ««<««« 

CHKR1  :  1 
CASE  IS  >  180 
CHKR1  :  0 

SATRAK  =  SATRAK  -  180 
END  SELECT 
SELECT  CASE  TGTHDG 
CASE  IS  <  0 
CLS 

LOCATE  12,  1 

PRINT  "*««*««««««*«*««*««****«*«**t<*««**««»*««» 

PRINT  "ERROR  ERROR  ERROR  -  DATA  IS  OUT 

PRINT  " - Re-enter  satellite  track - 

PRINT  " - Re-enter  target  heading - 

PRINT  ''«<**«*«*«*»«*«>*****««««*«**«*<***«**««*« 
CHKR2  :  1 
CASE  IS  <:  180 
CHKR2  =  0 


OFF  RANGE 


*<«*(*«**( 


««««***«** 

OFF  RANGE 


«*«***«*« 

ERROR  ERROR  ERROR  -  DATA  IS  OUT  OFF  RANGE 

- Re-enter  satellite  track - 

- Re-enter  target  heading - 

tt*t*t****t*****t*t*%**t*xttt*tttt*tttttt*t*ttt***t 
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CASE  IS  >  360 
CIS 

LOCATE  12,  1 

PRINT  "««««*«««<««««««««<«»«««««<««««<«*«**«<******«*«»**«' 

PRINT  “ERROR  ERROR  ERROR  -  DATA  IS  OUT  OFF  RANGE 


PRINT  " - Re-enter  satellite  track- 

PRINT  " - Re-enter  target  heading- 


PRINT  «««««»«" 

CHKR2  =  1 
CASE  IS  >  180 
CHKR2  =  0 

TGTHOG  =  TGTHOG  -  180 
END  SELECT 

BETA  =  ABS(SATRAK  -  TGTHOG)  /  57.29577951* 

ENO  FUNCTION 


■ 


I 


I  : 


SUB  CALC 

N02ANG  =  0NO2ANG  /  57.29577951* 

RNG  =  ( RANGE (SATALT))  *  1000 

T  =  TRANSMIT 

S2  s  SPOTS2(RNG,  BMOVG) 

IF  TPAREA  >  S2  THEN 
CHI  =  S2  /  TPAREA 

PRINT  "CHI  IS  CURRENTLT  AT— . CHI 

ELSE  CHI  =  1 
ENO  IF 

TLCS  =  LCS(N02ANG) 

CCT  =  (PG  »  JOULES  »  TLCS  *  RXAREA  «  (T  '  2)  »  CHI) 

CCS  =  (BMOVG  '  2)  *  (RNG  *  4) 

PWX  =  CCT  /  (PW  «  CCB  »  9.869604404*) 

S  s  (CCT  «  ML)  /  (CCB  *  1.96058E-23) 

’  PRINT  "o******************************* »<**<««<**«***«*«*«*««*«* 

’  PRINT  "proc  gain - PG 

’  PRINT  "joules - ";  JOULES 

’  PRINT  "rxarea - RXAREA 

’  PRINT  "transmittance-";  T 

’  PRINT  "chi - ";  chi 

’  PRINT  "beam  divergence";  BMOVG 

’  PRINT  "RANGE - “;  RNG 

’  PRINT  "pulse  width - ";  PW 

’  PRINT  PWX,  S 

END  SUB 

-f-f +-►++-► -f-*- 

'Significant  portions  of  this  program  were  structured  after 
’  interfaces  designs  built  or  described  by  Steven  Spence. 

’  AFIT/GST  89J 


SUB  CHGMENU 

'This  subroutine  changes  the  input  values  displayed  in  the  menu 
SELECT  CASE  SLCTR 
CASE  1 
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LOCATE  1,  33 
PRINT  "xxxxxxxx" 
LOCATE  1,  33 
INPUT  SATALT 
SELECT  CASE  SATALT 


CASE  IS  <  120 
CLS 

PRINT  " - - - WARNING  WARNING  WARNING -  " 

PRINT  • -  THIS  IS  PRETTY  LOW  -  " 

PRINT 

CASE  IS  >  30000 
CLS 

PRINT  - WARNING  WARNING 

— WARNING - ■■ 

PRINT  '■ -  THIS  IS  VERY  HIGH  FOR  A  TRACKING  SATELLITE 


PRINT 
CASE  ELSE 
PRINT 

END  SELECT 
CASE  2 

LOCATE  2,  33 
PRINT  ‘xxxxxxxx’ 

LOCATE  2,  33 
INPUT  LASER 
SELECT  CASE  LASER 
CASE  1 

WL  =  10.59104  /  1000000 
PRINT 

CASE  2 

WL  =  11.1494  /  1000000 
PRINT 

CASE  3 

WL  :  2.06  /  1000000 
PRINT 

CASE  4 

WL  =  1.057  /  1000000 
PRINT 
CASE  ELSE 
CLS 

COLOR  0,  2 
LOCATE  12,  28 
PRINT  "INPUT  ERROR" 

COLOR  2,  0 
LOCATE  24.  9 

INPUT  "HIT  RETURN  TO  CONTINUE  -  THEN  FIX  YOUR  ERROR";  PR$ 
END  SELECT 

CASE  3 

LOCATE  3,  33 
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PRINT  "xxxxxxxx" 

LOCATE  3,  33 
INPUT  JOULES 
CASE  4 

LOCATE  4,  33 
PRINT  "xxxxxxxx" 

LOCATE  4,  33 
INPUT  PWI 

PW  =  PWI  /  1000000 
CASE  5 

LOCATE  5,  33 
PRINT  "xxxxxxxx" 

LOCATE  5,  33 
INPUT  SNdB 

SNR  :  10  "  (SNdB  /  10) 

CASE  6 

LOCATE  6,  33 
PRINT  "xxxxxxxx" 

LOCATE  6,  33 
INPUT  RXmlcro 

RXSENS  =  RXmicro  /  1000000 
CASE  7 

LOCATE  7,  33 
PRINT  "xxxxxxxx" 

LOCATE  7,  33 
INPUT  PGdB 

PG  =  10  '  (PGdB  /  10) 

CASE  8 

LOCATE  8,  33 
PRINT  "xxxxxxxx" 

LOCATE  8,  33 
INPUT  RXD 

RXAREA  =  (3.14159265#  *  (RXO  ^  2 ) )  /  4 
CASE  9 

LOCATE  9.  33 
PRINT  "xxxxxxxx" 

LOCATE  9,  33 
INPUT  BEAMOO 
SELECT  CASE  BEAMOD 
CASE  1 
CLS 

LOCATE  12.  12 

INPUT  "ENTER  BEAM  DIVERGENCE  (In  Radians)" 
CLS 
CASE  2 

BMOVG  :  1.08  *  WL  /  RXD 


c 


■ 

CASE  3 

BHDVG  =  1.3  «  WL  /  RXO 
CASE  ELSE 

CLS 

COLOR  0,  2 

.  LOCATE  12,  28 

PRINT  "INPUT  ERROR" 

COLOR  2,  0 
LOCATE  24.  9 

INPUT  "HIT  RETURN  TO  CONTINUE  -  THEN  FIX  YOUR  ERROR';  PR| 
END  SELECT 

GaInTX  =  12.56637062#  /  (BMOVG  *  2) 

CASE  10 

LOCATE  10,  33 
PRINT  "xxxxxxxx" 

LOCATE  10,  33 
INPUT  ATMOS 
SELECT  CASE  ATMOS 
CASE  1.  2 

CASE  3 

^  LoLtE  7,  1 

■  PRINT  ■■  ENTER  THE  ATTENUATION  TO  ADD  TO  THE  ALREADY  EXISTING" 

PRINT  '•  ATTENUATION  FROM  THE  MOLECULAR  ABSORPTION  (in  dB)" 

PRINT  "  i.e.  50  X  attenuation  from  an  aerosol  would  be  3dB" 

PRINT  "  entry  is  for  ONE-WAY  passage  through  the  atmosphere" 

PRINT  "  at  a  zenith  angle  of  ZERO  (straight  up  through  the 
— atmosphere)" 

P  PRINT  "  — an  exponential  curve  (trans  vs  range)  is  used  for  all" 

PRINT  "  zenith  angles  other  than  zero.  ' 

INPUT  ATTENdB 

ATTEN  =  10  '  (ATTENdB  /  10) 

CASE  ELSE 
CLS 

COLOR  0,  2 

—  LOCATE  12,  28 

PRINT  "INPUT  ERROR" 

COLOR  2,  0  . 

LOCATE  24,  9 

INPUT  "HIT  RETURN  TO  CONTINUE  -  THEN  FIX  YOUR  ERROR";  PR$ 
.  END  SELECT 

CASE  11 

LOCATE  11,  33 
PRINT  "xxxxxxxx" 

LOCATE  11,  33 
INPUT  LCSMODEL 
i  CLS 

SELECT  CASE  LCSMODEL 
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CASE  1 

LOCATE  12.  5 

INPUT  "WHAT  IS  THE  RADIUS  OF  THE  SPHERE  IN  METERS" ;  RADSPHR 
CASE  2 
LOCATE  9,  4 

PRINT  "THE  FLYING  8RICK  MODEL  IS  A  THREE  DIMENSIONAL  MODEL" 

PRINT  "OF  AN  OBJECT  SIMPLIFYING  THE  TARGET  TO  THREE  ORTHOGONAL" 
PRINT  "PLANES  LIKE  THREE  SIDES  OF  A  BRICK.  THEREFORE,  THE  THREE" 
PRINT  "FACES  OF  THE  SYMETRICAL  BRICK  REPRESENT  THE  PROJECTED 

AREA" 

PRINT  "OF  THE  TARGET  FROM  A  FRONTAL/REAR  VIEW,  FROM  A  TOP  VIEW 

AND" 

PRINT  "FROM  A  SIDE  VIEW  ...  essentially  the  target’s  silhouette 
INPUT  "ENTER  THE  FRONT  /  REAR  PROJECTED  AREA  (sf4r  meters)  RA1 
INPUT  "ENTER  THE  TOP  PROJECTED  AREA  (sqr  meters) RA2 
INPUT  "ENTER  THE  SIDE  PROJECTED  AREA  (sqr  meters)";  RA3 
CASE  ELSE 
CLS 

COLOR  0,  2 
LOCATE  12,  28 
PRINT  “INPUT  ERROR" 

COLOR  2,  0 
LOCATE  24,  9 

INPUT  "HIT  RETURN  TO  CONTINUE  -  THEN  FIX  YOUR  ERROR";  PR| 
END  SELECT 

CASE  12 

LOCATE  12,  33 
PRINT  "xxxxxxxx" 

LOCATE  12,  33 
INPUT  REFLECT 
R  =  REFLECT  /  100 

CASE  13 

LOCATE  13,  33 
PRINT  "  xxxxxxxx" 

LOCATE  13,  33 
INPUT  SATRAK 
XXXX  =  BETA 

IF  CHKR1  :  1  THEN 
COLOR  2,  0 
LOCATE  24,  9 

INPUT  "HIT  RETURN  TO  CONTINUE  -  THEN  FIX  YOUR  ERROR";  PR$ 

END  IF 

CASE  14 

LOCATE  14,  33 
PRINT  "xxxxxxxx" 

LOCATE  14.  33 
INPUT  TGTHDG 
XXXX  =  BETA 

IF  CHKR2  =  1  THEN 


c 


II 

COLOR  2.  0 

LOCATE  24,  9 

INPUT  “HIT  RETURN  TO  CONTINUE  -  THEN  FIX  YOUR  ERROR";  PR$ 

END  IF 

CASE  15 

i  •  LOCATE  15,  33 

PRINT  “xxxxxxxx" 

LOCATE  15,  33 
INPUT  OFFSET 
CASE  16 

LOCATE  16,  33 
PRINT  "xxxxxxxx" 

LOCATE  16,  33 
INPUT  OATASLCT 

END  SELECT 
MENU 
END  SUB 


■ 


FUNCTION  dBoi  (PWX) 

.  V  s  PWX  *  1000 
IF  V  <  IE-09  THEN 
dBm  :  0 
GOTO  SKIPPER 
END  IF 

dBm  =  10  *  (LOG(V)  /  LOG(10«)) 
— 10 


SKIPPER: 

END  FUNCTION 


’X  is  in  watts 


’convert  to  milliwatts 

’take  loglO  of  v  and  mutiply  by 


LoglO  function  from  Microsoft 
QuickBASIC  4.0  BASIC  Language 
Reference,  1987 


FUNCTION  LCS  (N02ANG) 

SELECT  CASE  LCSMOOEL 
CASE  1 

LCS  =  R  *  (RAOSPHR  *  2)  *  8.377580414  ’pi*(8/3) 

CASE  2 

VL0»  =  RANGE (SATALT) 

VLm  =  XI  /  VLO* 


Sig1«  =  4  ♦  R  *  A1  *  ((VL14) 
Z  :  ZORNG  *  COS(NOZANG) 

VL2«  =  Z  /  VL0« 


Slg2i  =  4  ♦  R  *  A2  *  ((VL21I) 
VL34  =  OFFSET  /  VLOi 
Sig3«  s  4  *  R  *  A3  *  ((YL3i) 
LCS  =  slgl»  +  sig2»  +  sigSU 


CASE  ELSE 

PRINT  '■ - *»*♦**»  ERROR 

END  SELECT 


2) 


2) 

2) 


IN  FUNCTION  LCS 


*t***ttt 


i 
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END  FUNCTION 


SUB  MENU 
TRY AGAIN: 

CIS 

LOCATE  1,  1 

PRINT  "1.  SATELLITE  ALTITUDE" 
LOCATE  2.  1 

PRINT  "2.  LASER  TYPE  " 

LOCATE  3,  1 

PRINT  "3.  PULSE  ENERGY" 

LOCATE  4,  1 

PRINT  “4,  PULSE  WIDTH" 

LOCATE  5,  1 

PRINT  “5.  REQUIRED  S/N" 

LOCATE  6,  1 

PRINT  "6.  MIN  DETECTABLE  SIGNAL" 
LOCATE  7,  1 

PRINT  "7.  PROCESSING  GAIN" 

LOCATE  8,  1 

PRINT  "8.  DIAMETER  OF  OPTICS" 
LOCATE  9,  1 

PRINT  "9.  BEAM  DIVERGENCE  MODEL" 
LOCATE  10,  1 

PRINT  "10.  ATMOSPHERE  MODEL" 
LOCATE  11,  1 

PRINT  "11.  TARGET  MODEL" 

LOCATE  12,  1 

PRINT  "12.  TARGET  REFLECTANCE" 
LOCATE  13,  1 

PRINT  "13.  SATELLITE  HEADING" 
LOCATE  14,  1 

PRINT  "14.  TARGET  HEADING" 

LOCATE  15,  1 

PRINT  "15.  OFFSET  DISTANCE" 

LOCATE  16,  1 

PRINT  "16.  DATA  OUTPUT  REQUEST" 

LOCATE  1,  33 
PRINT  SATALT 
LOCATE  2,  33 
PRINT  LASER 
LOCATE  3,  33 
PRINT  JOULES 
LOCATE  4,  33 
PRINT  PWI 
LOCATE  5,  33 
PRINT  SNdB 
LOCATE  6,  33 
PRINT  RXmlcro 
LOCATE  7,  33 


19 


I 

PRINT  PGdB 
LOCATE  8,  33 
PRINT  RXO 
LOCATE  9.  33 
PRINT  BEAHOO 
im  LOCATE  10.  33 

PRINT  ATMOS 
LOCATE  11,  33 
PRINT  LCSMOOEL 
LOCATE  12.  33 
PRINT  REFLECT 
LOCATE  13,  33 
PRINT  SATRAK 
LOCATE  14,  33 
PRINT  TGTHDG 
LOCATE  IS,  33 
PRINT  OFFSET 

n  LOCATE  18,  33 

PRINT  OATASLCT 
LOCATE  19,  1 
COLOR  0,  2 

PRINT  "YOU  MAY  CHANGE  ANY  OF  THE  ABOVE  VALUES  SIMPLY  BY  ENTERING  ITS  NUMBER  ' 
COLOR  2,  0 

■PRINT 

PRINT  "ENTER  0  TO  CALCULATE  AGAIN  WITH  THE  ABOVE  DATA  SET" 

PRINT  "ENTER  21  IF  YOU  HAVE  HAO  ENOUGH  AND  WISH  TO  RETURN  TO  DOS" 

INPUT  SLCTR 
SELECT  CASE  SLCTR 
CASE  0 

EXIT  SUB 

■  CASE  21 

EXIT  SUB 
CASE  1  TO  16 

CHGMENU 
CASE  ELSE 

GOTO  TRYAGAIN 

_  END  SELECT 

END  SUB 

FUNCTION  RANGE  (SAT ALT) 

VR1*  :  6371.23  +  .01  ’ASSUMES  TARGET  ALTITUDE  OF  10  METERS 

•  VR2*  =  6371.23  SAT  ALT 

VR3*  =  (2  *  VR1#)  *  COS(3. 14159265#  -  NOZANG) 

VR4#  =  (VR2#  '  2)  -  (VR1#  *  2) 

VR5#  =  SQR((VR3#  '  2)  +  (4  *  VR4#)) 

’  ZORNG  Is  the  Zero  Offset  RaNGe  (i.e.  the  range  when  OFFSET  =  0) 

ZORNG  =  (VR5#  +  VR3«)  /  2 
I  XI  =  ZORNG  *  SIN(NOZANG) 

I  Z  =  ZORNG  *  COS (NOZANG) 

RNGSQR#  =  OFFSET  '2+X1“2+2'2 
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RANGE  =  SOR(RNGSQR») 

VR6i  =  SQR(RNGSQRi  -  (Z  *  2)) 

VR7i  =  VR6i  /  Z 
ZANG  =  ATN(VR7i) 

END  FUNCTION 

FUNCTION  SPOTSZ  (BMOVG,  RNG) 

SPOTSZ  =  .78539816»  *  (RNG  ‘  2)  *  (BMDVG  *  2) 
END  FUNCTION 

FUNCTION  TPAREA 


This  function  computes  the  total  projected  area  of  the  brick  model 
The  total  Is  a  sum  of  the  parts  (top, side  and  front) 


IF  LCSMOOEL  =  2  THEN 
0  =  RANGE (SAT ALT) 

API  =  XI  *  A1 
AP2  =  Z  *  A2 
AP3  s  OFFSET  *  A3 
SUM  =  API  +  AP2  +  AP3 
TPAREA  =  SUM  /  0 
END  IF 

IF  LCSMODEL  =  1  THEN  ! 

TPAREA  =  (RADSPHR  '  2  »  3. 14159265*) 

END  IF 

END  FUNCTION 
FUNCTION  TRAKTIME 

i 

VR1*  =  6371.23  +  .01  ’ASSUMES  TARGET  ALTITUDE  OF  10  METERS 

VR2*  =  6371.23  +  SATALT 

VR3»  =  (2  «  VR1*)  *  COSO.  14159265*  -  NOZANG) 

VR4#  =  (VR2*  ■  2)  -  (VR1#  ~  2) 

VR5*  :  SQR((VR3*  '  2)  ^  (4  «  VR4*)) 

’  ZORNG  Is  the  Zero  Offset  RaNGe  (i.e.  the  range  when  OFFSET  =  0) 

ZORNG  =  (VR5*  +  VR3*)  /  2 


TREV  =  (.009952  *  (VR2*  '  1.5))  /  60 
OMEGA  :  6.28318531*  /  TREV 
VK1  =  ZORNG  /  VR2* 

VK2  =  SIN( 3. 14159265*  -  NOZANG) 

VK3  =  VK1  «  VK2 
’  Calculate  the  ARCSIN  of  VK3 

VK4  :  SQRd  -  (VK3  '2)) 

— Tables 

BHALF  =  ATN(VK3  /  VK4) 

TRAKTIME  =  (2  *  BHALF)  /  OMEGA 
END  FUNCTION 


REV  Time  in  minutes 


'This  relationship  is  from 
'CRC  Standard  Mathematical 

’21st  edition,  page  233 


108 


G 


I 


II 


m 
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FUNCTION  TRANSMIT  Computes  transmittance  of  the 

’  atmosphere 

SELECT  CASE  LASER 


CASE  1  ’coa 

SELECT  CASE  ATMOS 
CASE  1,  3 


C1# 

z 

-.26438815# 

C2# 

r 

4.57337687# 

C3# 

= 

-8.2907455# 

C4# 

6.42929668# 

C5# 

= 

-1.88496131# 

Cl# 

= 

-.01635833# 

CASE  2 

Cl# 

z 

-.31883222# 

C2# 

z 

4.27150562# 

C3# 

z 

-7.28951884# 

C4# 

z 

5.38389615# 

CS# 

z 

-1.51345106# 

Cl# 

z 

-.0098359# 

END  SELECT 
CASE  2 

SELECT  CASE  ATMOS 
CASE  1,  3 

C1«  =  1.9240306# 

C2«  :  -4.72236837# 

C3#  =  6.25641669# 

C4#  =  -4.24200153# 

C5#  =  1.15335803# 

Cl#  =  .59438551# 

CASE  2 

Cl#  =  3.57172183# 

C2#  =  -8.747844260000001# 
C3#  =  11.56313808# 

C4#  :  -7.80928598# 

C5#  =  2.11201757# 

Cl#  =  .23738812# 

END  SELECT 

CASE  3 

SELECT  CASE  ATMOS 


CASE  1, 

3 

Cl# 

z 

.76297032# 

C2# 

z 

-1.96013817# 

C3# 

z 

2.66230004# 

C4# 

z 

-1,82902376# 

C5# 

z 

.50041431# 

Cl# 

z 

.85239347# 

CASE  2 

Cl# 

z 

3.98595311# 

C2# 

z 

-9.59816461# 

C3# 

z 

12.56345419# 

C4# 

z 

-8.43270371# 

usin9  C12 

’Clear  atmosphere 
'Molecular  absorption  only 


’HAZE  with  5  km  visibilty 

C02  with  C13 

’Clear  atmosphere 
’Molecular  absorption  only 

’HAZE  with  5  km  visibilty 


’Holmium 


’Clear  atmosphere 
’Molecular  absorption  only 


’HAZE  with  5  km  visibilty 
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CSi  =  2.27119836* 

Cl*  =  .12015229* 

END  SELECT 
CASE  4 

SELECT  CASE  ATHOS 
CASE  1,  3 

Cl*  :  .8226326* 

C2*  =  -2.09733802* 

C3*  :  2.83434822* 

C4«  =  -1.93997169* 

C5*  =  .52917037* 

Cl*  :  .83868596* 

CASE  2 

Cl*  =  3.67542735* 

C2*  =  -6.6265573* 

C3*  =  7.09392008* 

C4*  :  -4.1192832* 

C5*  :  .99615146* 

Cl*  =  -.25207037* 

ENO  SELECT 
END  SELECT 

VT1*  =  COS(ZANG) 

Q1  =  Cl*  *  (VT1*) 

Q2  =  C2*  *  (VTt*  *  2) 

Q3  =  C3*  «  (VT1*  '  3) 

Q4  :  C4*  »  (VT1*  ‘  4) 

Q5  =  C5*  «  (VT1*  ‘  5) 

IF  ATMOS  =  3  THEN 

TRANS  =  1  /  ATTEN 
ALPHA  =  ((LOG(TRANS))  /  120) 

L  :  RANGE! 120) 

TRANSA  =  EXP(ALPHA  *  L) 

ELSE 

TRANSA  :  1 

ENO  IF 

TRANSMIT  =  (Cl*  +  Q1  +  Q2  +  Q3  +  04  +  05)  »  TRANSA 
’  PRINT  ZANG,  VT1* 

ENO  FUNCTION 


’Neodyniun 


'Clear  atmosphere 
'Molecular  absorption  only 


'HAZE  with  5  km  visipilty 
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APPENDIX  B 


ROSEHA  STONE  FOR  OATS  MODEL 


OATS  MODEL 
VARIABLE 

A1 

A2 

A3 

ALPHA 
API 
AP2 
AP3 
ATMOS 
ATTEN 
ATTENdB 
BEAMOO 
BETA 
8HALF 

BMDVG 
BOTTOM 
Cl 

C2 

C3 

C4 

C5 


VARIABLE 
IN  THESIS 

~"'T  ■■■ 

At 

A2 

1 

A3 


a 


na 

na 

3 


VARIABLE 

DESCRIPTION 


transformed  area  of  the  frontal  surface 
the  brick 

transformed  area  of  the  too  surface  of 
the  brick 

transformed  area  of  the  side  surface  of 
the  brick 

attenuation  coefficient  for  additional  1 
mechan i sms 

frontal  component  of  the  target’s  total 
projected  area 

top  component  of  the  target’s  total 
projected  area 

side  component  of  the  target’s  total 
projected  area 

number  that  defines  the  atmosphere  type 
to  be  applied  to  the  model 

additional  atmospheric  attenuation  at  a 
zenith  angle  of  0  used  in  atmosphere 

the  variable  ATTEN  in  dB 


of 


OSS 


#3 


na 

na 

na 

na 

na 

na 

na 

na 


number  that  defines  the  beam  divergence 
type  to  be  applied  in  the  model 
angle  between  the  satellite  ground  track 
direction  and  the  target  heading 
geocentric  angle  between  the  target 

position  and  the  satellite  at  the  point 
that  the  satellite  is  first  trackable 
beam  divergence  angle  of  the  transmitted 
laser  beam  (full  angle  in  radians) 
beginning  value  for  the  variable  that  is 
to  be  incremented  when  using  DATASLCT  3 
first  coefficient  of  the  fifth  order 
polynomial  that  defines  atmospheric 
attenuation  for  the  model 
second  coefficient  of  the  fifth  order 
polynomial  that  defines  atmospheric 
attenuation  for  the  model 
third  coefficient  of  the  fifth  order 
polynomial  that  defines  atmospheric 
attenuation  for  the  model 
fourth  coefficient  of  the  fifth  order 
polynomial  that  defines  atmospheric 
attenuation  for  the  model 
fifth  coefficient  of  the  fifth  order 
polynomial  that  defines  atmospheric 
attenuation  for  the  model 


Cl 

CC8 

CCT 

CHI 

COUNT 

D 

DATASLCT 

DCHOICE 

DELTA 

DNOZANG 

GamTX 

INCH 

JOULES 

L 

LASER 

LCSMODEL 

LEFT 

MID 

N02ANG 

OFFSET 

OMEGA 

PG 

PGdB 

PW 

PWI 

PWX 


na 

na 

na 


attenuation  when  zenith  angle  is  zero, 

used  as  the  constant  in  the  fifth  order 
polynomial  defining  attenuation 
calculation  factor  in  subroutine  CALC 
IS  denominator  of  range  equation 
calculation  factor  in  subroutine  CALC 
is  numerator  of  range  equation 
ratio  of  laser  spot  size  to  target’s 
projected  area 

the  count  on  a  counter  in  the  routine  to 
find  the  zenith  angle  that  meets  the 
pre-establ ised  criterion 
the  value  in  function  TPAREA  that  defines 
the  range  from  target  to  satellite 


number  designating  type  of  data  to  be 
output  when  the  model  is  run 
number  designating  the  variable  to 
increment  for  DATASLCT  =  3 
the  difference  between  required  S/N  and 
the  actual  system  S/N  OR 
the  difference  between  detectable 
signal  power  and  actual  power 
the  no  offset  zenith  angle  measured  in 
degrees 

the  gain  in  the  system  from  the  laser 
beam  directionality 

the  value  by  which  the  variable  of  study 
is  incremented  when  DATASLCT  =  3 
the  energy  in  a  single  laser  pulse 
measured  in  joules 

length  of  atmosphere  through  which  a  laser 
beam  travels  -used  in  function  transmit 
number  designating  type  of  laser  to  be  used 
in  the  model 

number  designating  type  of  target  to  oe 
tracked  1  for  sphere  2  for  prick 
left  bound  for  the  halving  algorithm  used 
to  find  required  zenith  angle 
mid  value  for  the  nalving  algorithm  used 
to  find  required  zenith  angle 
no  offset  zenith  angle  (in  radians) 
distance  from  the  target  to  the  satellite 
ground  track 

revolution  rate  of  the  satellite 
processing  gain  of  the  receiver  system 
processing  gain  in  dB 

pulse  width  of  the  laser  pulse  in  seconds 
pulse  width  of  the  laser  pulse  in  i.is 
available  power  for  the  receiver  to  process 
(after  processing  gain  is  added) 
reflectance  of  the  target  skin 
raw  frontal  area  of  the  target 
raw  top  area  of  the  target 


R 

RA1 

RA2 


RA3 

A  3 

RAOSPHR 

na 

REFLECT 

na 

RIGHT 

na 

RNG 

d  or 

RNGSQR 

na 

RXAREA 

na 

RXD 

0 

RXmicro 

na 

RXSENS 

na 

S 

S/N 

SATALT 

h 

SATRAK 

na 

SIG1 

c\ 

S1G2 

02 

SIG3 

0-3 

SLCTR 

na 

SNdB 

na 

SNR 

na 

SUM 

na 

S2 

na 

T 

c 

TQTHOG 

na 

TLCS 

o\ 

TOP 

na 

TRANS 

na 

TRANSA 

Ts 

TREV 

P 

VK1 

na 

VK2 

na 

VK3 

na 

VK4 

na 

VLO 

na 

VL1 

na 

VL2 

na 

raw  side  area  of  the  target 
radius  of  the  target  sphere  (in  meters) 
reflectance  of  the  target  skin  in  percent 
right  bound  of  halving  algorithm  used  to 
find  the  required  zenith  angle 
range  of  the  target  from  the  satellite 
( in  meters) 

RNG  squared  (intermediate  calculation) 
effective  receiving  area  of  optics 
(in  square  meters) 

diameter  of  receiving  optics  (meters) 

minimum  detectable  signal  in  M-watts 

RXmicro  in  watts 

actual  S/N  at  the  receiver 

satellite  altitude  above  the  earth  (km) 


direction 

of 

satel 1 1 te 

ground 

track 

component 

of 

total 

LCS 

due 

to 

the 

transformed  frontal 

area 

component 

of 

total 

LCS 

due 

to 

the 

top  area 

component 

of 

total 

LCS 

due 

to 

the 

transformed  side  area 


required  S/N  for  target  declaration  in  dB 
required  S/N  for  target  declaration 
sum  of  area  components  for  projected  area 
computations 

soot  size  of  laser  beam  at  the  target 
transmittance  of  the  atmosphere 
target’s  heading  oriented  from  true  north 
target’s  laser  cross  section 
upper  bound  of  incremental  data  computation 
routine  when  DATASLCT  =  3 
transmittance  factor  of  the  atmosphere 
when  using  additional  attenuation 
transmittance  factor  which  modifies 
the  atmospheric  attenuation  when 
additional  attenuation  is  requested 


time  of  a  si 

ngle  sate 

1 1 i te  revo 

lution 

intermediate 

calculat 

1  on 

var 

1  ab 

le 

i  n 

the 

TRAKTIME 

function 

intermediate 

calculat 

ion 

var 

i  ab 

le 

i  n 

the 

TRAKTIME 

function 

intermediate 

calcu 1  at 

ion 

var 

1  ab 

le 

1  n 

the 

TRAKTIME 

function 

intermediate 

calculat 

ion 

var 

lab 

le 

i  n 

the 

TRAKTIME 

function 

intermediate 

calculat 

ion 

var 

lab 

1  e 

1  n 

the 

LCS  funct 

ion 

intermediate 

calculat 

ion 

var 

1  ab 

le 

1  n 

the 

LCS  funct 

ion 

intermediate 

calculat 

1  on 

var 

1  ab 

1  e 

1  n 

the 

LCS  function 
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VL3 


na 


VR1 

VR2 

VR3 

VR4 

VR5 

VR6 

VR7 

VT1 

WL 

XI 

Z 

ZANG 

ZORNG 


na 

na 

na 

na 

na 

na 

na 

na 

X 

X 

z 

e 

dt* 


ntermediate  calculation 

variable 

1  n 

the 

LCS  function 
ntermediate  calculation 

variable 

1  n 

the 

RANGE  function 
ntermediate  calculation 

variable 

1  n 

the 

RANGE  function 
ntermediate  calculation 

variable 

1  n 

the 

RANGE  function 
ntermediate  calculation 

variable 

1  n 

the 

RANGE  function 
ntermediate  calculation 

variable 

1  n 

the 

RANGE  function 
ntermediate  calculation 

variable 

1  n 

the 

RANGE  function 
ntermediate  calculation 

variable 

i  n 

the 

RANGE  function 
ntermediate  calculation 

variable 

1  n 

the 

TRANSMIT  function 
wavelength  of  the  laser 
X  component  of  distance  from  satellite 
to  the  target 

z  component  of  distance  from  satellite 
to  the  target 

zenith  angle  from  target  to  the  satellite 
range  from  satellite  to  target  if  OFFSET 
is  set  to  zero 


E 

APPENDIX  C  _  _  INPOT  DATA  FROM  FASOTE 
I  TnTot  ’  DATA'"FbF  C12  "’wifH~"NO"'’ HAZE 

0  cos  0  r  ( 0 ) 


0.000000 

1 .000000 

0.5465100 

3.000000 

0.9986295 

0.5460600 

5.000000 

0.9961947 

0.5452600 

8.000000 

0.9902681 

0.5432900 

10.00000 

0.9848077 

0.5414600 

12.00000 

0.9781476 

0.5392200 

15.00000 

0.9659258 

0.5350400 

17.00000 

0.9563047 

0.5317000 

20.00000 

0.9396926 

0.5258300 

23.00000 

0.9205049 

0.5188700 

25.00000 

0.9063078 

0.5136000 

28.00000 

0.8829476 

0.5046700 

30.00000 

0.8660254 

0.4980200 

33.00000 

0.8386706 

0 . 4868900 

35.00000 

0.8191521 

0.4786700 

37.00000 

0.7986355 

0.4697600 

40.00000 

0 . 7660444 

0 . 4550200 

43.00000 

0.7313537 

0 . 4384800 

45.00000 

0.7071068 

0 . 4263700 

47.00000 

0.6819984 

0.4133300 

50.00000 

0.6427376 

0.3918700 

53.00000 

0.6018150 

0.3679500 

55.00000 

0.5735765 

0.3505200 

57.00000 

0.5446391 

0.3318100 

60.00000 

0.5000000 

0.3012100 

63.00000 

0.4539905 

0.2674100 

65.00000 

0.4226182 

0 . 2430500 

67.00000 

0.3907312 

0.2172600 

70.00000 

0.3420202 

0.  1 761600 

73.00000 

0.2923718 

0. 1329800 

75.00000 

0.2588191 

0. 1039500 

77.00000 

0.2249510 

7.5797997E-02 

80.00000 

0. 1736482 

3.8364001E-02 

SELECTED 

PORTIONS  OF  SAS  OUTPUT 

R-SQUARE 


0.999995 

PARAMETER  ESTIMATE 


INTERCEPT 

A 

ASQ 

ACU 

AQUAD 

AOUIN 


-0.01635833 
-0.26438815 
4.57337687 
-8.29074550 
6.42929668 
-1 .88496131 


I 
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INPUT  DATA  FOR  C12  WITH  HAZE 


e 

cos  S 

T  id) 

0,00000 

1 .000000 

0.5238700 

10.00000 

0.9848077 

0.5186900 

20.00000 

0.93<^6926 

0.5026800 

30.00000 

O.860O254 

0.4742700 

40.00000 

0.7660444 

0.4305700 

50.00000 

0.6427876 

0.3669200 

60.00000 

0.5000000 

0.2767900 

70.00000 

0,3420202 

0. 1 556800 

80,00000 

0. 1736482 

3.0096000E-02 

SELECTED  PORTIONS  OF  SA3  OUTPUT 

PARAMETER 

R-SQUARE 

0.999999 

ESTIMATE 

INTERCEPT 

-0.00983590 

A 

-0.31883222 

ASQ 

4.27150562 

ACU 

-7,28951884 

AQUAD 

5.38389615 

AQUIN 

-1 .51345106 

'A  FOR  C13 


0.000000 

1 .000000 

0.9637000 

3.000000 

0.9986295 

0.9636500 

5.000000 

0.9961947 

0.9635600 

8.000000 

0.9902681 

0.9633500 

10.00000 

0.9848077 

0.9631500 

12.00000 

0.9781476 

0.9629000 

15.00000 

0.9659258 

0.9624400 

17.00000 

0.9563047 

0.9620700 

20.00000 

0.9396926 

0.9614100 

23.00000 

0.9205049 

0.9606300 

25.00000 

0.9063078 

0.9600200 

28.00000 

0.8829476 

0.9589900 

30.00000 

0.8660254 

0.9582000 

33.00000 

0.8386706 

0.9568700 

35.00000 

0.8191521 

0.9558700 

37.00000 

0.7986355 

0.9547600 

40.00000 

0.7660444 

0.9528800 

43.00000 

0.7313537 

0.9507000 

45.00000 

0.7071068 

0.9490600 

47,00000 

0.6819984 

0.9472300 

50.00000 

0.6427876 

0.9441 100 

53.00000 

0.6018150 

0.9404200 

55.00000 

0.5735765 

0.9375900 

57.00000 

0.5446391 

0.9343900 

60.00000 

0.5000000 

0.9287500 

63.00000 

0.4539905 

0.9218300 

65.00000 

0.4226182 

0.9162900 

68.00000 

0,3746066 

0.9061000 

70.00000 

0.3420202 

0,8976600 

73.00000 

0.2923718 

0.8814100 

75.00000 

0.2588191 

0.8671800 

77.00000 

0.2249510 

0.8488900 

SELECTED 

PORTIONS  OF  SAS 

R-SQUARE 

0.999971 


PARAMETER 


ESTIMATE 


INTERCEPT 

A 

A  SO 
ACU 
AQUAD 
AQUIN 


0.59438551 
1 .92403060 
-4.72236837 
6.25641659 
-4.2420153 
1 .15335803 


I 
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D 


i  INPUT  DATA  FOR  C13  WITH  RAZE 


e 

cos  e 

r  (0) 

0.00000 

1 .000000 

0.9269900 

10.00000 

0.9848077 

0.9259100 

20.00000 

0.9396926 

0.9224900 

30.00000 

0.8660254 

0.9161800 

40.00000 

0.7660444 

0.9057800 

50.00000 

0.6427876 

0.8887700 

60.00000 

0.5000000 

0.8593800 

70.00000 

0.3420202 

0.8014100 

80.00000 

0. 1 736482 

0.6476000 

I  I 


I  t 


SELECTED  PORTIONS  OF  SAS  OUTPUT 


R-SQUARE 

0.999996 


PARAMETER 

INTERCEPT 

A 

ASQ 

ACU 

AQUAD 

AQUIN 


ESTIMATE 

0.23738812 
3.57172183 
-8.74784426 
1 1 .56313808 
-7.80928598 
2. 1 1201757 


! 
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INPUT  DATA  FOR  HOLMIUM  WITH  NO  HAZE 


9  cos  6  t  (9) 


0.00000 

1 .000000 

0.9888800 

10.00000 

0.9848077 

0.9887100 

20.00000 

0.9396926 

0.9881 700 

30.00000 

0.8660254 

0.9871700 

40.00000 

0.7660444 

0.9855100 

50.00000 

0.6427876 

0.9827600 

60.00000 

0.5000000 

0.9779000 

70.00000 

0.3420202 

0.9679000 

80.00000 

0. 1 736482 

0.9381300 

SELECTED  PORTIONS  OF  SAS  OUTPUT 


PARAMETER 

R-SQUARE 

0.999992 

ESTIMATE 

INTERCEPT 

0.85239347 

A 

0.76297032 

ASQ 

-1 .96013817 

ACU 

2.66230004 

AQUAD 

-1 .82902376 

AQUIN 

0.50041431 
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INPUT  DATA  FOR  HOLMIUM  WITH  HAZE 


0 

COS  0 

T  (0) 

0.000000 

1 .000000 

0.9097400 

10.00000 

0.9848077 

0.9084100 

20.00000 

0.9396926 

0.9P42300 

30.00000 

0.8660254 

0.8965300 

40.00000 

0.7660444 

0.8838500 

50.00000 

0.6427876 

0.8631900 

60.00000 

0.5000000 

0.8277200 

70.00000 

0.3420202 

0.7586200 

80.00000 

0.1 736482 

0.5813500 

SELECTED  PORTIONS  OF  SAS  OUTPUT 


PARAMETER 


INTERCEPT 

A 

ASQ 

ACU 

AQUAD 

AQUIN 


R-SQUARE 

0.999997 

ESTIMATE 


0.12015229 
3.9859531 1 
-9.59816461 
12.56345419 
-8.43270371 
2.27119836 
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INPUT  DATA  FOR  NEODYMIUM  WITH  NO  HAZE 


6 

cos  d 

r  (6) 

0.000000 

1 .000000 

0.9874900 

10.00000 

0.9848077 

0.9873000 

20.00000 

0.9396926 

0.9867000 

30.00000 

0.8660254 

0.9855800 

40.00000 

0.7660444 

0.9837100 

50.00000 

0.6427876 

0.9806300 

60.00000 

0.5000000 

0.9751900 

70.00000 

0.3420202 

0.9640500 

8O.OC0OO 

0.1736482 

0.9314500 

SELECTED  PORTIONS  OF  SAS  OUTPUT 


PARAMETER 


INTERCEPT 

A 

A  SO 
ACU 
AQUAD 
AQUIN 


R-SQUARE 

0.999992 

ESTIMATE 


0.83868596 
0.82263260 
-2.09733802 
2.83434822 
-1 .93997169 
0.52917037 
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INPUT  DATA  FOR  NEODYMIUM  WITH  HAZE 


e 

cos  0 

T  (0) 

0.00000 

1 .000000 

0.7675500 

10.00000 

0.9848077 

•  0.7644300 

20.00000 

0.9396926 

0.7546400 

30.00000 

0.8660254 

0.7367900 

40.00000 

0.7660444 

0.7080200 

50.00000 

0.6427876 

0.6627100 

60.00000 

0.5000000 

0.5893700 

70.00000 

0.3420202 

0.4619700 

80.00000 

0. 1 736482 

0.2199000 

SELECTED  PORTIONS  OF  SAS  OUTPUT 

R-SOUARE 

1 .000000 


PARAMETER 

ESTIMATE 

INTERCEPT 

-0.25207037 

A 

3.67542735 

ASQ 

-6.62655730 

ACU 

7.09392008 

AQUAD 

-4.1 1928320 

AQUIN 

0.99615146 

<•*  > 


APPENDIX  D 


«  SAMPLE  FASCODE  INPUT  DATA  FILES  * 
C12~at  40  tegrcM  zenith  angltt 


C03/nohaz*/KZSUmiXtlTER/40I»S 

; 

1  ^ 

HI-1  F4-1 

CN-1  AE-1 

01-1  SC-0 

FI-0  PL-0 

TS-0  AT-1 

MG-0  LS-1 

0 

0  0 

944.194 

944.194 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

j 

0.000 

0.000 

1 

3  3 

0  1 

1  26 

0 

0.000 

120.000 

944 . 194 

0.010 

0.000 

40.000 

0.000 

0.000 

0 

•* 

0.000 

0.000 

0.000 

0.000 

0.000 

0  0 

0  3 

0  0 

0.000 

0.000 

0.000 

0.000 

0.000 

1 

C13-zt  40  Oeoreez  zenith  angle 

BZ613C03/M0KAZE/KZDLATWZNTER/4  OOEGREES 

HI-1  F4-1 

CN-1  AE-1 

01-1  SC-0 

FI-0  PI/-0 

TS-0  AT-1 

MS-0  LS-1 

0 

0  0 

696.909 

896.909 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

1  t. 

3  3 

0  1 

1  28 

0 

0.000 

120.000 

896.909 

i 

f 

0.010 

0.000 

40.000 

0.000 

0.000 

0 

1 

0.000 

0.000 

0.000 

0.000 

0.000 

0  0 

0  3 

0  0 

0.000 

0.000 

0.000 

0.000 

0.000 

HOIMIUH-at 

40  degrees  zenith  angle 

1  I 

HOLKIOK  IN  YLr/NOHAZE/^!lZajlT«nNTER/40DECREES 

i 

r  ^ 

HI-1  F4-1 

CN-1  AE-1 

EM-1  SC-0 

PI-0  PL-0 

TS-0  AT-1 

MG-0  LS-1 

0 

0  0  < 

4854.369 

4654.369 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

r  ■ 

0.000 

0.000 

3  3 

0  1 

1  28 

0 

0.000 

120.000 

4854.369 

0.010 

0.000 

40.000 

0.000 

0.000 

0 

1 

0.000 

0.000 

0.000 

0.000 

0.000 

i  . 

0  0 

0  3 

0  0 

0.000 

0.000 

0.000 

0.000 

0.000 

j 

F  ■ 

NEODYNIUM-at  40  degrees  zenith  angle 

1 

— NBODYNIUM/NO  HAZE/MIDLATITODE  H1NTER/40DEG— 

Hl-1  F4-1 

CN-1  AE-1 

EM-1  SC-0 

FI-0  PL-0 

TS-0  AT-1 

HS-0  LS-1 

0 

0  0 

9460.736 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

1 

3  3 

0  0 

1  28 

0 

0.000 

120.000 

9460.738 

1 

0.010 

0.000 

40.000 

0.000 

0.000 

0 

0.000 

0.000 

0.000 

0.000 

0.000 

r 

0  0 

0  3 

0  0 

0.000 

0.000 

0.000 

0.000 

0.000 
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